LIMIT OPERATORS FOR CIRCULAR ENSEMBLES 

KENNETH MAPLES, JOSEPH NAJNUDEL, AND ASHKAN NIKEGHBALI 

Abstract. It is known that a unitary matrix can be decomposed into a product of 
reflections, one for each dimension, and the Haar measure on the unitary group pushes 
forward to independent uniform measures on the reflections. We consider the sequence 
of unitary matrices given by successive products of random reflections. 

In this coupling, we show that powers of the sequence of matrices converge in a 
suitable sense to a flow of operators which acts on a random vector space. The vector 

CO ' space has an explicit description as a subspace of the space of sequences of complex 

numbers. The eigenvalues of the matrices converge almost surely to the eigenvalues of 
the flow, which are distributed in law according to a sine-kernel point process. The 
eigenvectors of the matrices converge almost surely to vectors which are distributed in 
law as Gaussian random fields on a countable set. 

.^ , This flow gives the first example of a random operator with a spectrum distributed 

according to a sine-kernel point process which is naturally constructed from finite di- 
mensional random matrix theory. 
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Notation 

If f G C" is a vector, then we write v[m\ for the image of v under the canonical 
projection map C" — ?> C™ onto the first m standard basis vectors and we write {v)k for 
its k-th coordinate. 

We write 0{n) for the orthogonal group of dimension n, i.e. the group of invertible 
operators on M" which preserve the standard real inner product. We write U{n) for 
the unitary group of dimension n which preserves the standard complex inner product. 
For any vector space V , real or complex, we write GLiV) for the group of invertible 
transformations. We always write 1 for the identity operator in every space. 

We write U = U{1) for the unit circle in C, i.e. those complex numbers with modulus 1. 

Calligraphic characters denote a-algebras, i.e. A, B, C, etc. If A and B are a-algebras 
on a common set, then AV B denotes the smallest u- algebra containing both A and B. 

We also write a\/ b, for a,b > 0, to mean max(a, b) and a A 6 to mean min(a, b). 

We employ asymptotic notation for inequalities where precise constants are not impor- 
tant. In particular, we write X = 0{Y) to mean that there exists a constant C > such 
that |X| < CY. All constants are assumed to be absolute unless indicated otherwise 
by an appropriate subscript; e.g. we write fn{x) = On{gn{x)) to mean that there are 
constants C„, one for each value of n, such that |/„(a;)| < C„5'„(a:) for all x. We also use 
(modified) Vinogradov notation, where X <Y means X = 0{Y), for convenience. 

If t is a real number, we write [tj its integer part. 

If /J is a Hilbert space with scalar product (.,.), and ii F G H, then F-^ = {x G 
H; {x,y) = Vy G F}. li H is a complex Hilbert space, then we will always use the 
scalar product which is linear in the first variable and conjugate linear in the second, 
i.e. {ax, by) = ab{x,y). 

1. Introduction 

It has been observed that, for many models of random matrices, the eigenvalues have 
a limiting short-scale behavior when the dimension goes to infinity, which depends on 
the global symmetries of the model, but not on its detailed features. For example, the 
Gaussian Orthogonal Ensemble (GOE), for which the matrices are real symmetric, with 
independent gaussian entries on and above the diagonal, corresponds to a limiting short- 
scale behavior for the eigenvalues, which is also obtained for several other models of 
random real symmetric matrices. Similarly, the limiting spectral behavior of a large class 
of random hermitian and unitary ensembles, including the Gaussian Unitary Ensemble 
(GUE, with independent, complex gaussians above the diagonal), and the Circular Uni- 
tary Ensemble (CUE, corresponding to the Haar measure on the unitary group of a given 
dimension), involves a remarkable random point process, called determinantal sine-kernel 
process. This process can be defined as follows: it is a point process for which the fc-point 
correlation function is given by 



Pfc(xi, . . . , Xk) = det ( \ ^ "__ '" 

TXyXp Xqj 



l<:P,q<k 



From an observation of Montgomery, it has been conjectured that the limiting short- 
scale behavior of the imaginary parts of the zeros of the Riemann zeta function is also 
described by a determinantal sine-kernel process. This similar behavior supports the 
conjecture of Hilbert and Polya, who suggested that the non-trivial zeros of the Riemann 
zeta functions should be interpreted as the spectrum of an operator | + iH with H an 
unbounded Hermitian operator. 
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In order to understand the kind of randomness which would be involved in such an 
operator, it is natural to try to construct a random version Hq of it, for which the spectrum 
has the conjectured limiting behavior of the zeros of the Riemann zeta function, i.e. is a 
determinantal sine-kernel point process. Since the spectrum of Hq should also correspond 
to the limiting behavior of many ensembles of hermitian and unitary matrices, it is natural 
to expect that these ensembles can, in one way and another, be related to Hq. 

Instead of looking directly for Hq, one can also directly seek the flow of linear operators 
(t/^)aeM := (e*"^°)aeR generated by exponentiation. This point of view is both consistent 
with what would be a possible interpretation in quantum mechanics {Uq playing the role 
of the operator of evolution at time a, whereas Hq corresponds to the Hamiltonian) , and 
with the number theoretic point of view: if the Hilbert-Polya operator H would exist, 
the Chebyshev function ipQ, which is defined as 

v;^ 1 I log p, X = p™ for some prime power p^ 

~-^ I U, otherwise 

where the summation is through all powers of primes bounded by x, would formally 
satisfy 

7/'o(e^) = /" (e^ - e^/2 Tr(e*^^)) dy + 0(1), 

J —oo 

which suggests an important role played by the conjectural flow [e^^'^^y^^ of unitary 
operators. To see this, recall the von Mangoldt formula (see e.g. [T5]): 

^o(a;)=a;-5^- + 0(l), 
p 
where the summation is over all zeros p of the Riemann zeta function. Now taking x = e^ 

in the above, writing e^ — ^ — as the integral of its derivative and then writing the 

last sum as the trace of the operator yields the formal identity above. 

It should be mentioned that the problem of the existence of the operator Hq is also 
hinted at in the work by Katz and Sarnak [6J. For instance, the zeta function of algebraic 
curves are related to the unitary group or some other compact groups (e.g. the symplectic 
group) and there the question of coupling all different dimensions of the unitary group 
together in a consistent way in order to prove strong limit theorems (i.e. almost sure 
convergence) and having an infinite dimension space and operator sitting above is raised. 

The main goal of the present paper is the construction of a flow of random operators 
(l^")Q,giR5 whose spectrum, in a sense which can be made precise, is a determinantal 
sine-kernel process, and which is directly related to the Circular Unitary Ensemble. A 
similar, but simpler, construction has been made in [IT], in which we consider permutation 
matrices instead of unitary matrices. The construction which is made in the present paper 
uses the recent construction of virtual isometries given in ^ and needs several further 
steps. 

First, the space on which (V^")Q,gM acts must be infinite-dimensional. In order to relate 
this space with the Circular Unitary Ensemble, we will construct a coupling between the 
matrix models of each dimension. That is, we define a sequence ('U„)n>i of random unitary 
matrices, in such a way that for all n > 1, each matrix m„ G f/(n) is Haar distributed on 
the unitary group. Of course, there are many ways to couple the random variables m„ to 
each other: for example, by taking all the matrices to be independent. However, in order 
to have sufficient consistency to construct limiting objects, we need to be more subtle. 

In fact, we will couple (m„)„>i in such a way that almost surely, this sequence of random 
matrices is a virtual isometry. The notion of virtual isometry has been introduced in one 
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of our former articles [2] , generalizing both the notion of virtual permutation studied by 
Kerov, Olshanski, Vershik [S], and the previous notion of virtual unitary group introduced 
by Neretin JT2\. By definition, a virtual isometry is a sequence («„)„>! of unitary matrices 
Un of dimension n, such that for all n, Un is the matrix u such that for Un+i fixed the 
rank of 



Q 1 / ^n+1 

is minimal (which, as we will see, is always or 1). From this definition, one can deduce 
that for all n > 1, Un determines completely ui,U2, ■ ■ ■ ,Un-i, and from these matrices, 
one directly obtains a decomposition of m„ as a product of complex reflections. Moreover, 
the Haar measures in different dimensions are compatible with respect to the notion of 
virtual isometrics, i.e. it is possible to construct a probability distribution on the space of 
virtual isometrics in such a way that for all n > 1, the n-dimensional component follows 
the Haar measure on U{n). Therefore, the Circular Unitary Ensemble can be coupled in 
all dimensions by considering a random virtual isometry following a suitable probability 
distribution. 

In [2], it is shown that if {u„)n>i follows this distribution, then for all k, the kth 
positive (resp. negative) eigenangle of m„, multiplied by n/27r (i.e. the inverse of the 
average spacing between eigenangles for any matrix in U{n)), converges almost surely 
to a random variable yk (resp. yi-k)- The random set {yk)kez is a determinantal sine- 
kernel process, and for each k, the convergence holds with a rate dominated by some 
negative power of n. In the present paper, we improve our estimate of this rate, and more 
importantly, we prove that almost sure convergence not only holds for the eigenangles 
of Un, but also for the components of the corresponding eigenvectors. More precisely, we 
show that, for all k,i > 1, the ^th component of the eigenvector of u„ associated to the 
kth positive (resp. negative) eigenangle converges almost surely to a non-zero limit tk/ 
(resp. ti-k/) when n goes to infinity, if the norm of eigenvector is taken equal to ^/n and 
if the phases are suitably chosen. Moreover, the variables {tk,i)k£Z,e>i are iid complex 
gaussians. 

Knowing the joint convergence of the renormalized eigenvalues and the corresponding 
eigenvectors, it is natural to expect that, in a sense which has to be made precise, the 
limiting behavior of {un)n>i when n goes to infinity can be described by an operator 
whose eigenvectors are the sequences (tk/)e>i, k E Z and for which the corresponding 
eigenvalues are {yk)k&- 

The precise statement of our main result can in fact be described as follows: 

Theorem 1.1. Almost surely there exists a random vector subspace T of C°° and a flow 
of linear maps {V°')aeR on J-", such that V""^^ = V^V^ , and satisfying the following 
properties, for any sequence w = {wi)e>i in T : 

(1) Foranyfixedi> I, the ith component of then-dimensional vector of Un{wi, ... ,Wn) 
tends to the ith component ofV^^w) when n goes to infinity. 

(2) Ifwy^O, the L^ distance between ulC {wi, . . . ,Wn) and {{V"w)i, . . . , {V"w)n) is 
negligible with respect to the norm of {wi, . . . , Wn). 

Moreover, the eigenvectors of the flow, i.e. the sequences w E J^ such that there exists 
A G M for which V°'w = e^*'^-^"^ for all a E M., are exactly the sequences proportional to 
{tk,e)i>i for some k E'L. The corresponding value of X is equal to yk. 

Intuitively, the operator V"" can be viewed as a limit, for sequences in the space J-", of 
the iteration Un \ when n goes to infinity. The flow {V°')aeR can be compared with the 
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flow constructed in pG[j for permutation matrices. Tlie space T is random and explicitly 
defined in terms of the virtual isometry (Mn)n>i- 

Some natural questions which can then be asked are the following: 

(1) Is it possible to replace the space T by another space with a simpler description? 

(2) Is it possible to define a version of the flow (V"")^^^ of operators, which can be 
naturally related to ensembles of unitary or hermitian matrices, which are different 
from the CUE? 

An answer to this last question would give a more generic version of (V^")agR, which may 
enlighten in a new, more geometric way, the properties of universality enjoyed by the 
sine-kernel process in random matrix theory. 

Such a generalization seems to be very difficult to construct. In particular, we do not 
know how one could couple the GUE in such a way that the renormalized eigenvalues 
converge almost surely. Note that the successive minors of an infinite GUE matrix cannot 
converge almost surely to a limiting distribution. 

It would also be interesting to relate the operators we construct to the Brownian 
carousel construction, which was introduced by Valko and Virag in [12], in order to 
generalize the sine-kernel process to the setting of /3-ensembles for any /3 G (0, oo) (instead 
of just /3 = 2). 

The paper is devoted to providing the details of the construction of the fiow of operators 
given in Theorem 11.11 Our construction is both geometric and probabilistic and this is 
refiected in our arguments which involve for instance coupling techniques or martingale 
convergence theorems. For the convenience of the reader, we will first briefiy outline 
in the next section, without proof, the construction of the fiow of operators (\^")aeK in 
Theorem 11.11 

2. Overview of the construction 

In Section E] we define the set of complex reflections of a finite dimensional complex 
vector space. These refiections are unitary maps which are linear combinations of the 
identity and a map of rank one. It turns out that a unitary matrix u G U{n) can be 
written as a product 

(1) u = rn---ri 

where each r^, 1 < /c < n is a complex refiection. In this decomposition, each reflection r^ 
fixes every standard basis vector Cj with k < j < n, so that the partial products r^ ■ ■ ■ ri 
are naturally unitary maps on U{k). This decomposition has already been used in pLj to 
provide a new approach to the study of the characteristic polynomial of random unitary 
matrices (in particular in relation with the moments conjecture of Keating and Snaith 
[7]) and it was also exploited in [3] in the study of the circular Jacobi ensemble. 

Decomposition ([T]) was also used in [2] to define virtual isometries {un)n>i as families 
of unitary matrices, m„ G U{n), so that «„ = r„M„_i for every n > 1 for some family of 
refiections {r„,)n>i- The push-forward of Haar measure on U{n) decomposes in a simple 
way: the distribution of each r^ is independent, and each rk{ek) takes values uniformly 
on the unit sphere in C'^. The Haar measure on U{n) for n finite can then be naturally 
extended to a uniform measure on the set of virtual isometries. 

The next step is to find an explicit relation between the eigenvalues and eigenvectors 

of Un+i and u„ for each n > 1. Let us write Xi = e*^i , . . . , An = e*^" for the 
eigenvalues of m„, almost surely distinct, with the ordering < 6'|"' < ■ ■ ■ < 6n < 27r. 
Let fl"' G C" denote a unit length representative of the eigenspace for A^" . Then if we 
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expand r„+i(e„+i) in the basis of eigenvectors 

n 

i=i 
then the eigenvalues of m„+i are precisely the zeros of the rational equation 

" \W I -I 1 2 

El (")|2 '^3 ^ \^ ~ ^'^\ 1 - 

3=1 ^3 '^ 

and the eigenvectors of Un+i are given by the n + 1 equations 

n (n) 

^ .(n+1) ^ -^ ^i ^W , ^n- ^ 

for 1 < /c < n + 1; here C^ G M^ is a constant so that /^" has unit length. 

In Section O we make the following key observation: the sequence of eigenvalues A^" , 

with \ <k <n and ra > 1, is independent of the argument of the coefficients /Xj"^ , with 
1 < J < 'T- and n > 1. Therefore, we can consider the sequence of eigenvalues of the 
virtual isometry and prove that it converges almost surely, and then, conditioning on the 
eigenvalues of every matrix in the virtual isometry, consider the sequence of eigenvectors 
and show that they also converge in a suitable sense. 

The first part of this plan is carried out in Section [71 where it is shown that the 
eigenangles converge to a limiting process {yk)kez in the sense that 

^^"^ =yk + 0,((1 + A;2)n-H^) 

almost surely for every e > 0. 

Next, in Section [8] we condition on the eigenangles of the entire sequence of matrices. 
We show that for each fixed k there is a renormalization factor Dj!^ so that for each i > 1 
the sequence {D]^ f]^\ei/) is a martingale which converges in L^ and almost surely to a 
limiting value g^/- For each A; > 1, we call the sequence {gk,i)i>i € C°° the eigenvector 
of the isometry. 

We require stronger information on the convergence of the eigenvectors in order to 
make a suitable definition. In Section [10] we consider just those sequences in C°° which 
are spanned by the eigenvectors of the isometry. Then we show that we can define a 
"fiow" {U")aeR on this vector space by considering sequences Un t[n], where (t^)£>i is 
a linear combination of eigenvectors. In fact, for each eigenvector g^ and £ > 1 we show 
that {un gk[n],ei) converges almost surely for every a to e^^^^^^^gk/- Furthermore, we 
show for all 5 > that uniformly for a in compact sets, 

so that the fiow converges strongly (in this sense) in the limit. 

We are now ready to define our random vector space and the associated fiow. We define 
our random space F to be those sequences {wi)i>i, along with sequences {{y'^w)()i>i, so 
that Un w converges weakly to V°'w and so that 



\u. 



i""Jw[n]-(y"«;)[n]|| =0{n^~ 



for some 5 > and choice of the implied constant. This latter inequality ensures that 
we do not introduce additional "exotic" eigenvalues in the limit. We then define the fiow 
V^ : w H-> V°'w from the definition of J-". 
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3. Reflections and the definition of virtual isometries 

In this section, we recall the definition and elementary properties of refiections and 
the way they are used to construct virtual isometries. The reader can refer to |2] for 
more details but we wanted to give an independent presentation here to make the paper 
self-contained and also to fix notation and conventions. The idea behind is to provide 
a way to generate the Haar measure with the help of "simple" and independent unitary 
transformations. The idea is not new and a for instance a general algorithm is given by 

3.1. Refiections over M. We will begin by briefiy recalling the definition of refiections 
over R since this is the one which is used if one wants to carry our construction to the 
orthogonal group. We also wish to state them here in order to understand why real 
refiections or Householder transformations would not be suitable if the ground field is the 
field of complex numbers. 

Definition 3.1. Let H denote a real vector space of dimension n. A refiection of H 
(sometimes called a Householder transformation or an elementary refiector) is 

defined as an orthogonal transformation of H which fixes each element of a hyperplane 
(i.e. linear subspace of codimension 1). 

Let r be a refiection. Let us assume that r 7^ 1 so that ker(l — r) = Kisa, hyperplane 
and dimim(l — r) = 1. 

Since r G 0(]R) we must have det r G {±1}. Clearly 1 is an eigenvalue with multiplicity 
n — 1. Let A 7^ 1 denote the other eigenvalue. As detr is the product of the eigenvalues 
we conclude A = — 1 and det r = — 1. Consequently there exists a vector a E H such that 
r{a) = —a and note that r is a map of order 2 (namely r^ = 1). We will write r^ for the 
refiection which maps a 1— )■ —a. 

It is easy to find a formula for r^: First note that Ma © (Ma)-*" = H. There exists 
G H* (where H* is the space of linear forms on H) such that x — ra{x) = (f){x)a for 
every x E H and ker = ker(l — r^,) = K. Since ker((-, a)) = K, there exists a non-zero 
A G C* such that = A(-,a). Therefore x — ra{x) = X{x,a)a. Now ra{a) = —a so we 
have A = 7-^, hence 

{a,aj ' 

(2) raix) =x- 2p^a. 

[a, a) 

This is the equation of a real orthogonal refiection that sends a to —a and hence a 
householder transformation is parametrized by a vector a. 

Finally, we note that ii m,e E H are two distinct vectors with ||?7i|| = ||e|| = 1, there 
exists a unique refiection which maps m on e, namely r^-e- But it is easy to see that 
given two vectors of norm 1 in a complex Hilbert space, then there does not necessarily 
exist a Householder transformation ([2]) which maps one on the other. Hence we need to 
introduce another type of refiection. 

3.2. Refiections over C Over the complex numbers the theory is different since the 
reflections we consider are not necessarily of flnite order and they are parametrized by 
one vector and a phase (that is, a complex number of modulus 1). 

We now assume we are given a Hilbert space H with dim if = n. For u,v E H, we use 
the standard inner product 

n 

{u,v) := y^^UkVk- 

k=l 

If F C if , F a subset of H, we write F-^ := {x G ii | {u, x) = for all u G F}. 
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If F and G are subspaces of H ^ we write H = F ©^ G to indicate H = F + G = 
{x + y\xGF,y& G} and (x, y) = for all x G F and y E G. It is then easy to check 
that F = G^ and (G^)^ = G and dimG + dimG^ = dimi/ = n. 

Let U{H) denote the set of unitary operators, i.e. linear bijections u : H ^ H which 
preserve the inner product: (ux, uy) = {x, y) for every x,y E H. Let 1 denote the identity 
map. We first recall the elementary and well known result: 



Lemma 3.2. If u E U{H) then im(l — u) = ker(l — u 



,± 



Definition 3.3. Let u denote a linear transform of a complex vector space H of finite 
dimension. We call u a complex reflection if it is the identity or if it is unitary and 
rank(l — u) = 1. We will just write reflection if the base field is clear. 

Remark 3.4. Usually in the literature, a linear transformation g G GL{H) is a reflection 
if the order of g is finite and rank(l — u) = 1. For our applications we do not require that 
g have finite order. 

As in the real case, we can compute a formula for a general complex reflection. 

Proposition 3.5. Suppose that r is a reflection of the space H and that the vector a E H 
spans im(l — r). Then there exists a G U such that for every x E H, 

r[x) = X — (1 — a)- ;-a. 

(a, a) 

Proof. We have by construction and Lemma 13.21 the equivalences im(l — -u) = Ca and 
ker(l — u) = (Ca)-*". Let G H* denote the linear form defined by u{x) = x — 0(x)a. 
Clearly 0(x) = if and only if (1 — u)x = 0, so ker0 = ker(l — u). But the linear 
form (■, a) also vanishes on ker(l — u) = (Ca)-*-, so we must have some A G C* such that 
(j) = X{-,a). Hence m(x) = x — A(x, a)a. 

To determine A, we note that u{a) = aa for some a G C, \a\ = 1, so we must have 
A = -7^ as required. D 

Definition 3.6. For non-zero a E H and a G C, |a| = 1, we define the reflection raQ,(x) 
by 

ra,a{^) = X - (1 - a) ' a. 

[a, a) 

Note that r^ „ has eigenvalue 1 with multiplicity n—1 and eigenvalue a with multiplicity 
1. The following facts are easy to check and we omit the proofs. 

Proposition 3.7. For any non-zero a E H and a, (3 eV we have the following. 

(1) ra,ara,0 = Ta^ap- 

(2) For every g E U{H), gVa^^g* = rga,a- 

(3) For every non-zero A G C, rxa,a = i^a,a- 

(4) r-^a = ra,a- 

Proposition 3.8. Let a,b E H be non-zero vectors and a, (3 E V. Then the reflections 
Ta^a o,nd Th^fB commutc if and only if Ca = Cb or (a, b) = 0. 

Proof. This follows immediately by writing 

ra,arbA^) = X - 1 - «H-^a - 1 - /? ^if + ^ ~ " ^ ~ /3 \//, « 

{a, a) {b,b) {a,a}{b,b} 

which shows that the reflections commute if and only if 

{b,a){x,b)a = {a,b){x,a)b. D 
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Now we note that given two distinct vectors e,m & H oi unit length, there exists a 



1 — (m,e) 
/ DUL,ii tiidt / \^o; — /;t, vviiiL,ii lo i m~e,a vviicic u: — 

r is given by the equation 



unique complex reflection r such that r(e) = m, which is rm~ea where a = — ^^^ Such 

^ ' ' 1 — (m,e) 



(x, m-e) 

r[x) = X [m — e). 

1 — {m, e) 

3.3. Virtual isometries and projections of unitary matrices. We now explain how 
to construct an infinite dimensional structure which contains in a natural way each finite 
dimensional unitary group. It is crucial to our construction that the unitary matrices 
between different dimensions be closely linked. This will allow us to give a meaningful 
definition to almost sure convergence of random matrices as the dimension grows to 
infinity. 

Proposition 3.9. Let H he a complex Hilbert space, E a finite dimensional suhspace of 
H and F a suhspace of E. Then for any unitary operator u acting on H which fixes 
every vector in E-^ , there exists a unique unitary operator 'Ke,f{u) on H which satisfies 
the following two conditions. 

(1) 'Ke,f{u) fixes every vector in F-^ D E-^ . 

(2) The image of H under u — tte,f{u) is contained in the image of F^ under u — 1. 
Moreover if G is a suhspace of F, then ttf,g ° '^e,f{u) is a well-defined unitary operator 
on H and is equal to 't^e,g{'^)- 

Proof. First we prove uniqueness. Let a; G -F fl (m — 1)(-F-'-). There exists y G F-^ such 
that X = u{y) — y. Since x E F and y G F-^ we have 

My)f = \\yf + \\xf 

by the Pythagorean theorem. Since u is unitary ||'u(t/)|| = ||y|| so we would require 
||x|| = 0, hence F n {u - 1)(F^) = {0}. 

Now if Vi and V2 are two unitary operators satisfying the properties of tte,f{u), then: 

(1) Vi and V2 fix globally F since they fix F-^. 

(2) Vi — V2 vanishes on F-^ by construction. 

(3) The range of vi — V2 is included in (m — 1){F-^) since the range of each of f i — -u 
and u — V2 are. 

We conclude from (2) and (3) that the image oivi—V2 is contained in F\^{u—1){F-^) = {0} 
and so the operators must agree. 

Next we prove the tower property of 'n'E,F{u), assuming that the operator exists. 

Let G C F C E C H. We write v = nE,F{u) and w = 71e,g{v). These operators are 
well-defined and satisfy 

(1) V fixes every vector in F-^ 

(2) {u-v){H)C{u-l){F^) 

(3) w fixes every vector in G-*- 

(4) {v-w){H) C (v-l){G^). 
This yields the elementary calculation 

{u - w){H) C span{(M - v){H), {v - w){H)} 

Cspan{(tz-l)(F^),(t;-l)(G^)} 

C span{(M - 1)(F^), {u - 1){G^), {u - v){G^)} 

C span{(M - 1){G^), {u - v){H)} 

C{u-1){G^). 
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Since w fixes each vector of G-^, it is equal to 7!'e,g{'^)- 

Now we prove tlie existence of the operator by induction. It is sufficient to prove the 
existence of tie^f in the particular case where E = span{F, e}, where e is a unit vector 
orthogonal to F. In this case, if m is a unitary operator fixing each vector of E-^, then 
the operator v = tte,f{u) can be constructed explicitly as follows. 

If u{e) = e then on taking v = u which fixes spanjii^-'-, e}, hence it fixes F-^ and 
{u-v){H) = (m-1)(F^). 

If u{e) 7^ e then for all x E H we define 

{u{x),e-u{e)) ^ 

v{x) := u{x) r~r\ — \^ — \^ " ""l^JJ = r o u 

1- {u{e),e} 

where r, as indicated, is the unique refiection mapping u{e) i-)- e. Hence w is a unitary 
transformation. 

Now let X G -E"*" so that u{x) = x and e — ^(e) G E since E is globally fixed by u. Here 
(m(x), e — u{e)) = {x,e — u{e)) = so v{x) = x. Moreover by construction v{e) = e, so 
consequently v fixes each vector in F^. 

Finally, for all x E H,we have u{x)—v{x) = '~f{e—u{e)) for some 7 G C, so {u—v){H) C 
{u — 1)(F-'-), hence v satisfies the conditions of tte,f{u). D 

Remark 3.10. It follows from Proposition 13.71 that if r := (r)~^ is the refiection mapping 
e onto u{e) then 

u = rTiE,F{u). 

Similarly one can easily prove that u = 7Te,f ° f'' where r' is the unique refiection such 
that r'{u~^{e)) = e. 

The existence of the projection map described above suggests how to define "virtual 
isometries", the infinite dimensional objects alluded to above. Indeed, let H = £^(C) 
and {ei)£>i be the canonical Hilbert basis of H. Then for all ra > 1 the space of unitary 
operators fixing each element of Vn := span(ei, . . . , e^)-*- can be canonically identified 
with the unitary group U{n). 

By identification, for n > m > 1 the projection 7rv„,v,n gives a map from U{n) to [/{m), 
more simply noted TTn,m, so that, for n > m > p > 1, 

Definition 3.11. A virtual isometry is a sequence {un)n>i of unitary matrices such 
that for all n > 1, m„ G U{n) and 7r„_|_i^„(u„+i) = m„. The space of virtual isometries will 
be denoted f/°°. 

Remark 3.12. U°° does not appear to have a natural group structure. 

Remark 3.13. In this definition we made a particular choice of vector spaces lying in £^(C). 
Although it appears to be necessary to make a choice for the probabilistic arguments later, 
the ideas in this construction apply also to other Hilbert spaces with a sequence of basis 
vectors. 

Proposition 3.14. Let (x„)„>i he a sequence of vectors with Xn G C" lying on the unit 
sphere (i.e. ||a;„||^2(c") = ^)- Then there exists a unique virtual isometry (m„)„>i such 
that Un{en) = Xn for every n > 1. In particular, Un is given by 

where for each I < j < n, rj = 1 ifxj = Cj and otherwise is the unique refiection mapping 

Cj I— )■ Xj. 
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Proof. This follows directly from the proof of Proposition 13.91 and the remarks which 
follow it. n 

Remark 3.15. In the particular case where every x„ = et(n) for some t{n) G {1, . . . ,n}, 
then (m„)„>i is the sequence of matrices associated to a virtual permutation [11]. The se- 
quence of permutations {<Jn)n>i is constructed by the so-called Chinese restaurant process 
for all n > 1, 



C"n — Tn,t{n)Tn-l,t{n-l) " " ' ''"1,1 

where Tkj = 1 ii j = k and otherwise is the transposition (j, k). 

Remark 3.16. The above proposition shows in particular the fact that any unitary matrix 
of U{n) can be expanded into a product of reflections. 

Proposition 13.141 shows in particular that U°° is non-empty. Moreover, our construction 
will allow us to construct measures on U°°. Indeed, it is natural to look for the analogue 
of Haar measure on ?7°°. First, we recall the following correspondence from |^ (this can 
also be obtained as a consequence of the results above and the subgroup algorithm of 
Diaconis and Shahshahani [1]). 

Proposition 3.17 ([2J). Let (x„)„>i be a random sequence of vectors with x„ G C" and 
let {un)n>i be the unique virtual isometry such that M„(e„) = x„ for all n > 1. Then 
for each n > 1 the matrix Un is distributed according to Haar measure if and only if 
Xi,...,Xn are independent and for each I < j < n, Xj is distributed according to the 
uniform measure on the unit sphere in . 

As a consequence, we deduce the compatibihty between Haar measure on U{n), n > 1, 
and the projections 7r„ ^ for n > m > 1. 

Proposition 3.18. For all n > m > 1 , the push-forward of the Haar measure on U{n) 
under 'Kn,ra is the Haar measure on U{m). 

Remark 3.19. One has to check that 7r„m is measurable as is done in [2] 

The compatibility property of Proposition 13.181 allows us to define the Haar measure 

on U°°. 

Proposition 3.20. Let U denote the a-algebra on U°° generated by the sets 

W{k,B):={{Un)n>l\UkeB} 

where k > 1 and B C U{k) is a Borel set. Let {fin)n>i be a family of probability measures, 
fin defined on {U{n),B{U{n))), where B{U{n)) is the Borel set on U{n), and such that 
the push-forward of fin+i by 7rn+i,n is equal to fin- Then there exists a unique probability 
measure /i on {U°°,U) such that push-forward of fi by the nth coordinate map is equal to 
Un for all n > 1. 

Proof. This is the Kolmogorov extension theorem applied to the family (/in)n>i- See 

H- " □ 

Now combining Proposition 13. 171 and Proposition 13.201 we obtain the following. 

Proposition 3.21 ([2]). There exists a unique probability measure fi^°°^ on the space 
{U°^,U) such that its push-forward by the coordinate maps is equal to the Haar measure 
on the corresponding unitary group. In particular, let {Xn)n>i be a random sequence of 
vectors with each Xn G C" almost surely on the unit sphere, and let {un)n>i be the unique 
virtual isometry such that M„(en) = x„. Then the distribution of {un)n>i is equal to iJ.^°°^ 
if and only if the x„ are independent random variables and for all n > 1, Xn is uniformly 
distributed on the unit sphere. 
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4. Spectral analysis of the virtual isometries 

It is classical that if m„ is a random unitary matrix following the Haar measure on 
U{n), then the distribution of the eigenangles of m„, multiplied by n/27r, converges in law 
to a determinantal sine-kernel process. 

In ^ it was proven that the eigenangles of a virtual isometry, taken according to Haar 
measure and renormalizing the eigenangles by the dimension n, converge almost surely 
to a point process with this determinantal distribution. Precisely, we have the following. 

Proposition 4.1. Let {un)n>i be a random virtual isometry following the Haar measure, 
and let 

■ ■■< e'^^l < e^^l < e^^^ < o < ^i"^ <9P <... 

be the eigenangles of Un- Then almost surely, for all m ^ Z, there exists ym such that 

when n goes to infinity, e > being some universal constant. Moreover, the process 
{ym)m&z is a determinantal sine-kernel process. 

To construct our limiting flow of operators we need to understand the behavior of the 
eigenvectors of m„ as n goes to infinity. Our method will give a formula for the spectrum 
and eigenvectors of m„+i in terms of the spectrum and eigenvectors of m„. 

We assume throughout that for each n > 1, the n eigenvalues of m„ are distinct; this 
holds almost surely for virtual isometries constructed according to the Haar measure in 
Proposition I3.21[ 

Clearly each matrix m„ has its eigenvalues supported on the unit circle. We write 
A^" , A2 , • • • , An for the eigenvalues of m„, ordered in such a way that if A^"^ = e*^* 
then almost surely under Haar measure, 

< ^S"^ < ■ ■ ■ < ^i") < 27r. 

(These inequalities define what is sometimes called a Weyl chamber of the group). More- 
over, the eigenangles enjoy a property of periodicity: for all k E Z, 0j^_ln = ^i" + 2vr. 

As all the eigenvalues are distinct, each eigenvalue corresponds to a 1 dimensional 
eigenspace. We can therefore write /}" , ..., /« for the family of unit length eigenvectors 
of Un, which are well-defined up to a complex phase. 

Let Xn = Un{en) and let r„ denote the unique reflection on C" mapping e„ to Xn- 
Therefore, we have Un+i = t„+i o («„ © 1). It is natural to decompose x„+i into the 
basis given by i(/}" ),..., 6(/„), e„+i, where t : C" — )■ C"'~^^ is the inclusion which maps 
(xi, . . . , Xn) to (xi, . . . , Xn, 0). Identifying fl and i(/fc ), we then have 

n 

Xn+1 = 2_^ ^^k fk + ^n^n+i 
fc=l 

for some yU^ (^ < k < n) and z/„ such that |/ii P + " " ' + l/^n P + Wn\'^ = 1- 

The following result gives the spectral decomposition of Un+i in function of the decom- 
position of Un and Xn+l'- 



Theorem 4.2 (Spectral decomposition). Conditionally on the event that the coefficients 
Hi , . . . , Hn are all different from zero and that the n eigenvalues of Un are all distinct 
(which holds almost surely under the uniform measure on U°° ), then almost surely the 
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eigenvalues of Un+i are the unique roots of the rational equation 



A,- i — V, 



|2 



on the unit circle. Furthermore, they interlace between 1 and the eigenvalues of u^, 

< ^!"+'^ < ^i") < ^5"+'^ < • ■ ■ < ^i") < ^i;V^ < 27r. 
anc? ^/le eigenvectors satisfy the relation 

n (n) 

/L(n+l)a An+l) _ Sr^ f^j An) ^n - ^ 

y''k ) Jk ~ A^ Un) .{n+l)J3 "^-, . (n+1) ^^"+1 

j=l ^j - ^k -L - A^ 



where 



, (") 1 2 



(n+l) _ \^ \^^j I I ^n ~ 1 1 

(n) \ (n+l) 12 H _ \(n+l)<2 



A — t ^ v"-7 \ V'-i-i; 2 h X*- 

j=l l^j - ^fc I l-L - ^fc 

is i/ie unique positive real which makes fj^ a unit vector. 

Proof. Let / be an eigenvector of Un+i with corresponding eigenvalue z. Then we have 

n 

i=i 

where ai,...,a„,6 are (as yet unknown) complex numbers, not all zero. Our goal is to 
write these coefficients in terms of x„+i and the eigenvalues of «„. 
We have 

zf = Un+lf 



y^ ttjUn+lfj + bUn+lCn+i 



i=i 

n 

i=i 
We recall that for all t G C"+\ rn+i(t) is given by 

r„+i(t) = t + r{Xn+i - e„+i) 

SO that 

zf = y^ ajK /, + 7 r(a;„+i - e„+i) + 6a;„+i. 

^ J \^ J (e„+i,a;„+i -e„+i) ' 

Now we decompose 

n 

E{n) r{n) 
f^k Jk ' ^n^n+l 
fc=l 
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and 



-/=E«iM"Mr + 



(n) I An) . ^^. 



(n) 



X 



n+l ~ ^n+l) 



bXn+1 



i=i 



E'-iM-vr + IE'"^ 

i=i 



(n) /^^ 



(n) 



/^^ - 1 



l^^n+l ~ 6n+l, 



+ 6x 



n+l- 



cin) 



fin) 



Because /{ , ..., fn , e„+i is a basis for C"+ , we deduce the system of n + 1 equations 



zttj = aj 



(n) 



Af+^«->5:.,A«=i^+vr 



Z^r, 



for j = 1, . . . ,n and 



z6 = 6 + (z/„ — 1) 2_] o.e^. 



[n) ^h 



in) 



e=i 



Z^r, - 1 



6(z/„ 



,W 



,W 



For z ^ {Aj^" , . . . , Xn , 1}, let us define the hnear transform Q 
matrix representation on the basis /{" , . . . , /« , e„+i is 



y+l _^ £n+l ^j^Qgg 



where 



Q = / + wv\ 



W 



,(") 



A^ — ^ 
\ 1-2 / 



and 



r„) M^"' 



,(«)Mi"' 



Then, the above system can be written 

Q/ = 0. 

Clearly, rankQ G {n, n + 1}. If it has full rank then / = 0, but we assume a priori 
that z is an eigenvalue for Un+i and so has a non-trivial eigenspace. Thus we must have 
rank Q = n and 

The right hand side can only vanish if / is proportional tow,sof = aw for some complex 
constant a E C\ {0} and v^w = — 1. In particular. 



n x(n)|,,(")|2 
>r^ A,- l/i. 



J lO 



An) 



+ 



Wn-M 
1-Z 



I - l^n 



3=1 ^i 
as required. 

Conversely, if 2; ^ {^i\ • • • ; ^n , 1} satisfies this equation, then 

Qw = w + wiy^w) = w + w{—\) = 0, 

which implies that w is an eigenvector of Un+i for the eigenvalue z. 

Let us now show that the eigenvalues z ^ {A^" , . . . , An , 1} of Un+i strictly interlace 
between 1 and the eigenvalues of m„: since Un+i has at most n+l eigenvalues, this implies 
that \i , . . . , An , 1 are not eigenvalues of m„+i. 



LIMIT OPERATORS FOR CIRCULAR ENSEMBLES 15 

Define the rational function ^iS*^— T-CUloojby 

Note tliat $ vanislies precisely on the eigenvalues of Un+i which are different from 
A-[ , . . . , An , 1. Recalling that |/^i P + • • ■ + |/^n T + knP = I5 '^^ can rearrange the 
expression of $ to the equivalent form 




Hence, $ takes values only in iM. U {00}, since for all z ^ z' ^ 5^, (2; + z')/{z — z') 
is purely imaginary (the triangle {—z\z^z') has a right angle at z). Note that for z G 
{A^" , . . . , An , 1}, a unique term of the sum defining $ is infinite, since by assumption, 
/i^"" , . . . , ^n -I ^ — i^n are nonzero and X{ \ . . . , An , 1 are distinct: <l>(2;) = 00. 

Next, we consider 1 1— )■ $(e**) in a short interval (6':" —5, ^:" +5). Then, for t = 6^ +u 
in this interval, 

4i f^ = T^— = 2m-i + 1 

while the other terms in $(e**) are uniformly bounded as 5 — )■ 0; likewise for the interval 
(—5, 5). In particular, $ — )■ ioc as m — )■ from the right and $ — )■ — icxD as m — )■ from 
the left. We therefore conclude, as $ is continuous, that on each interval of the partition 

(0, ^i"^) U (^!"\ ^^"^) U • ■ ■ U (0^"), 27r) 

of the unit circle into n+ 1 parts, 1 1— )■ $(e**) must assume every value on the line iM, and 
in particular must have at least one root. But we know that $ only has n + 1 roots on 
the circle so there must be exactly one root in each part of the partition, which proves 
the interlacing property. 

It remains to check the expression of the eigenvectors {f^ )i<k<n+i given in the 
theorem, but this expression is immediately deduced from the expression of the vector w 
involved in the operator Q defined above, and the fact that \\f^ || = 1. 

D 

5. A FILTRATION ADAPTED TO THE VIRTUAL ISOMETRY 

In our proof of convergence of eigenvalues and eigenvectors, we will use some martingale 
arguments. That is why in this section, we introduce a filtration related to our random 
virtual isometry. For ra > 1, we define the a-algebra An = <j{a'^ \ ^ "^ rn < n^l < j < 
m} and its limit A = V^^iAn- 

Lemma 5.1. For all n > 1, the a-algebra An is equal, up to completion, to the a-algebra 
generated by Ui the variables |/U^ | and Um for 1 <m <n — 1 and 1 < j < m. 



Proof. Under Haar measure, the assumptions of Theorem 14.21 are almost surely satisfied 
for all n > 1. Hence, almost surely, the eigenvalues of Un+i are the roots of the equation 

n \{") 1^ 1 2 
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This equation depends only on |/ij |, A^ for j = 1,...,?t,, and z/„. By induction, we 

deduce that A^" is a measurable function of ui, {fij }i<j<m,i<m<n and {i'm}i<m<n- 

Conversely, the above equation with z specialized to A^" , . . . , X^+i can be written 
in the form 

Rv = w. 
Here w is a column vector of Is, v is the column vector with entries 

y^ = I 1^1 I . . . l/^n I 1 — 1/ 



and i? is an (n + 1) x (n + 1) matrix with entries 

A(") 

3 ^ 

for 1 < j < n and 



^_^(n+l)' 



If we write R = RS, where S is the diagonal matrix with entries 5*^^ = {fij Xj )~^ 
(1 < j < n) and Sn+i,n+i = ^n, — ^, then we see that the rows of R are the representations 
of the eigenvectors /{" , • • • , /n+i in the basis /}" , . . . , /« , Cn+i up to constants, and 
therefore orthogonal. We conclude that R, and thus R, are invertible, so Rv = w has 
a unique solution, which can be written in terms of the eigenvalues X^ for 1 < j < m 

and m < n + 1. Thus we conclude that |/ij" | and z/„ are measurable functions of A^"^ for 
1 < j < m and m < n + 1 as was to be shown. D 

For 1 < j < n, we define the phase (f)j by fij = 0^ l/ij"^ |, and the a- algebras 
Bn = Ay (r{(f)f'^ I 1 < m < n - 1, 1 < j < m} and i3 = V^^^S^. 

Lemma 5.2. The a-algebra Bn is equal, up to completion, to the a-algebra generated by 
A and the eigenvectors / for I < j < m and 1 < m < n. 

Proof. Again by Theorem 14. 2^ we can write the eigenvectors of m„+i in function of fXj = 

^("■Ji^Wj (1 < j < n), Vni A, (1 < j < m, m G {n,n + 1}), and the eigenvectors 



'J lO \ \- — J — -Ji -n, -'J 

■'J 



of Un- Clearly, 1/^^ |, Vn, and A^ for m G {n,n + 1} are ^-measurable and cpf" is 
,B„_|_ i-measurable . 
Conversely, we write 

An) I (n) I 
/|'/,("+l)U/2 f(n+l) An). _ Vj l/^i I 
\y''k ) Jk ^Jj I (n) ,{n+l) 

^- ~ ^k 

to find each <\y^ as a function of the other variables. D 

6. A PRIORI ESTIMATES FOR UNITARY MATRICES 

In this section, we give some bounds on the probability of different events under the 
Haar measure on the unitary group. Let us fix e > 0, and let us define the following 
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events: 

Eq = {9^^^ ^ 0} n {\/n > 1, z/„ ^ 0} n {Vn > 1, 1 < A; < n, /ij^"^ ^ 0} 

El = {ziriQ > l,Vn > no, jz/^l < n~2+'=j. 

E2 = {3^0 > 1, Vn > no, 1 < A; < n, |/i["^| < n"5+^} 

^3 = {3^0 > 1, Vn > no, k > l,n-l'' < 4+i - 0^ < n~'^'}. 

We then let E := Eq n Ei n E2 Ci E3. We will show that E holds almost surely for each 
e > 0, so that we have a priori information on the distribution of the eigenvalues of the 
random virtual isometry. It will turn out that these estimates imply strong quantitative 
bounds on the change in eigenvalues of successive unitary matrices. 

Remark 6.1. In [2j an analogous event was defined to avoid small values of the dimension 
n where the behavior of the eigenvalues can be more erratic. This event will serve the 
same purpose, although we have chosen the exponents more carefully to sharpen our 
results. 



We begin by showing that for any fixed basis of C", the coefficients of a uniform random 
vector on the unit sphere are almost surely 0(n,~2+'^) for any e > 0. 



Lemma 6.2. Suppose fi,...,f„ G C" is an orthonormal basis and x G C", ||a;|| = 1 is 
chosen uniformly from the unit sphere. Then if we write x = xii^i + ■ ■ ■ + XnVn, we have 
the large deviation bound 



xj\'> 5) = 0(exp(-(5n/6)), 



for all 6 > and j = 1, 



Proof. Let us assume S < 1, since the probability is zero otherwise. Let 0i,...,0n be 
iid uniform random phases in S^ and let ei, ...,€„ be independent standard exponential 
random variable. Then, it is well-known that x has the same distribution as the random 
vector y = yiVi + ■ ■ ■ + ynVn where 



Vj = <P3 



ej 



ei H VCr. 



Now 



n 



so 



P(ei + ■ ■ ■ + e„ < -) < e"/2E[exp(-ei e„)] < e"/22-« < e-"/^ 

P(|xjf >5)< P(ej- > — ) + 0(exp(-n/6)) 

= 0(exp(-(5n/6)), 

since 5 G (0,1). D 

From this estimate, we deduce the following bound on the coordinates of the eigenvec- 
tors ff^. 

Lemma 6.3. Let e > 0. Then, almost surely, we have 



sup |(/f),e,)p = 0(n-i+^) 

l<j,(.<n 



and 



f{«) 



i(/r'e.)n^ 



0{n- „ 



sup E 
where the implied constant may depend on e and {um)m>i- 



i+e^ 
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Proof. Consider the vector / for each fixed j and n. Let .^ G C, |^| = 1 be an 
independent random phase. Then by the invariance by conjugation of Haar measure, the 

'j 



vector $,/)"' is distributed in law uniformly on the unit sphere. One deduces that for all 



¥{\{&\e,)\' > n-i+^) = 0(exp(-n76)). 



'J 
Using Borel-Cantelli lemma gives the first result. Moreover 



Jo 





oo 

'o 

oo 




We deduce: 

P(E[|(/f\e,)r/^|^] >n^-t) <nt~%[E[|(/f\e,)r/^|^] 

= n7-%[|(/f),e,)|«/^] = 0(n-^). 
By Borel-Cantelli lemma, for all but finitely many n > 1, 1 < j,i < n, 

E[|(/f,e,)rA|^]<n^-|. 
By the Holder inequality applied to the conditional expectation, for e sufficiently small. 

Another consequence of Lemma 16.21 is the following: 
Proposition 6.4. The events Eq, Ei, E2 all hold almost surely. 
Proof. We apply Lemma 16.21 to the decomposition 

n 

E,,in)An) 



which gives 

so, in particular. 



P(l/^i"^r > n-^+') = 0(exp(-n76)) 



n>l l<k<n 



Therefore, by the Borel-Cantelli lemma, almost surely only a finite number of the events 
{|/^fc T > n~^~^^} hold simultaneously. A similar argument controls the coefficients z/„. 

n 

Before we can control E^, we require some estimates on the eigenvalues of a Haar 
unitary random matrix. Recall that if m„ is distributed according to Haar measure, then 

(n) 

one can define, for 1 < p < n, the p-point correlation function pp of the eigenangles, as 
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follows: for any bounded, measurable function from MP to R, 



E 



E <^(^r.---'^S 



l(n) ^(n)^ 



P 

[0,27r)P 



pj,"^ (ti, . . . , tp)(f){ti, . . . , tp)rfti . . . dtp. 



Moreover, if the kernel K is defined by 

sin(rit/2) 



Kit) :- 



2tt sin(t/2) 
then the p-point correlation function can be given by 

p(")(ti,...,t„) = det(ir(t,-t,))^^^^^. 

Let us first show that the gaps between eigenvalues cannot be asymptotically much 
larger than average. 

Lemma 6.5. Let I C [0, 27r) be Lebesgue measurable. Then 

\F\ 
F{all of the eigenvalues of Un are in I) < exp( n). 

2tc 

Proof. We recall the Andreiev-Heine identity [14J, which says that 

P(all of the eigenvalues of u„ are in /) = det M 
where M' is an n x n matrix with entries 



Ml, = jeMKj-k)t)^ 



for j,k between 1 and n. Note that the matrix (exp(i(j — k)t)^ _ is hermitian and 

positive, M^ is also; likewise M^" . By the orthogonality of characters M^ + M^" = 1 so 
M^ , M^" have the same eigenvectors and, if we write r,- for 1, . . . ,n for the eigenvalues 
of M^" , then 1 — Tj are the eigenvalues of M^ . The eigenvalues of each matrix must lie 
in the interval [0, 1], as otherwise one of the eigenvalues of the other matrix would be 
negative. Now, 



n 



det M^ = TT(1 - Tj) < exp(- V r.) = exp(- Tr M^') = exp( n) 

i=i i=i 

as was to be shown. D 

Note that the previous lemma applies to all measurable subsets, although we will only 
need to apply it to intervals. 

Next we control the gaps between eigenvalues from below. 

Lemma 6.6. Suppose ti, ...,tp C / lie in an interval of length \I\ = S < 1/n. Then we 
have the estimate 

p(")(ti,...,g = Op(<52^-v^-2). 

Proof. We have 

p^;\t,,...,t^) = det{Kit,-t,))l^^^. 

The Taylor series for the sine function shows that for |t| < 1/n, 

Kit) = ^(l-^^{n'-l)t^ + 0{nH-)). 
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Thus, we have: 

pt'Htu • • • , ^n) = ^ det (1 - ^(n^ - l)(t,, - t,r + 0{n\U - t,r))l^. 

Let A denote the p x p matrix in the last display, let 1 denote the column vector of all 
ones and let Wj denote the column vector whose ith entry is 

K), = 1 - A, = ^(n^ - l)(t, - t,.)2 + Oin\U - t,f). 
Then by multilinearity and the inclusion-exclusion principle, 

det A = y^ (— l)''^' det (t"! ■■■ Up) , where w,- = < ■'' 

,^ 1, otherwise 

Clearly each term is zero if more than one of the columns is equal to 1 , so we get 

p 
det A = {-iy~^ Yl ^^* ^J + (-l)''det M 
i=i 
where M is the matrix with columns Wi, ...,Wp and Mj is M with the jth column re- 
placed with 1. Then in the expansion of each determinant we can bound each term by 
Op{{n^6'^y~^), and the conclusion follows. D 

Proposition 6.7. The event E^ holds almost surely. 

Proof. Fix n > 1. The probability that two adjacent eigenvalues of Un differ by at least 
2S is bounded above by the probability that one of the parts of the partition 

(0, 6)U{6,2S)U---U{ [2716-^ \ 6, [2716-^ + Ij 6) 

contains no eigenvalue. This, by Lemma 16. 5^ is bounded by 

[27i6-^ + IJ exp{-6n/27i). 

Now we let S = n~^~^'^ and apply the Borel-Cantelli lemma to show that at most a finite 
number of the Un have gaps larger than n~^~^'^. 

Next, we see by Lemma 16.61 that the probability that two adjacent eigenvalues of Un 
differ by at most S < 1/n is bounded by 

p'i'\ti,t2) dtidt2 = 0{n^6^) 

\tl~t2\<S 

which, when we specialize 6 = n~3~^^ is 0{n~^~'^)] summing over n and applying the 
Borel-Cantelli lemma shows that these events occur at must a finite number of times as 
well. D 

7. Convergence of the eigenangles 

Lemma 7.1. Let e > 0. Then, almost surely under the Haar measure on U°° , for n>l 
and < k < n^^^ , we have 

/d("+1)| ('^)|2 



and for n>l and —rt}!'^ < ^ < 0, 

^(n+l)| (n) 1 2 



^(") _ ^("+1) 



0((l + |A;|)n-3- 
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Remark 7.2. The implied constant in the O(-) notation depends on {um)m>i and e: in 
particular, it is a random variable. However, for given {um)m>i and e, it does not depend 
on k and n. 

Proof. By symmetry of the situation, we can assume /c > 0. Moreover, let us fix e G 
(0, 0.01): we now suppose that the event E occurs for this value of e, which is true almost 
surely. Recall from Theorem 14.21 that 



\(«)|,,(n)|2 ,, 

A- l/i- I |i — Z/„ 



|2 



j=i ^j - ^k ^~ ^k 



l^n.. 



By extending n-periodically the notation X"' , fi"' , / to any j G Z, we can write 



\("-)|,,(")|2 i-i 12 

V^ Aj- \Hj I |1 — Z/„| _ 

j6j Aj- - Afc -L - ^fc 



where J is the random set of n consecutive integers, such that 9^ — vr < 9^ < 9]^ ' +tt. 
Iterating the lower bound on the distance between adjacent eigenvalues, given by the 
definition of the event i^s, we get, for j G J\{k — 1, k}, 

\ef-et^'\>\k-3\n~i-\ 

and then 

\\(n) , (n+l)| ^ I , .1 -^-e 

l^i -K I ^ \k-3\n « \ 

smce |t^] c^A: I ^ ^• 

Likewise, we have by Eg, 1 - Ai"+^^ = 0(A;n-i+^), and by £"2, |/if^P = 0{n-^+'^'), 
which gives, for j G J\{k — 1, A;}, 

— ^ n 3^ 

^(n)_^(n+l) -\k-3\ 

Summing for j in J\{k — 1, A;}, which is included in the interval [fc — 1 — n, fc + n], gives 

jeJ\{fc-i,fc} ^3 ^k 

Now, subtracting this equation from the product of ([3]) by 1 — A^"" , and bounding 
i^n = 0(n~2+'=) (by the property Ei) gives us the resulting equation 

Next we estimate the first two terms in terms of the eigenangles. We find 

and 

Af) - Ai"+^) = z(ef ) - 4"+^))Af ) + O((0f ) - 9^^^y) 

ioT j = k — 1, k. Collecting terms and using the trivial bounds gives 



3 
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^(n+l), (n)|2 

]("+l)| ,,(") 12 






From Theorem 14.2^ the eigenvalues of Un and u„+i interlace, so for n sufficiently large 
the real part of the first term is positive and the real part of the second term is negative. 
The real part of the right hand side tends to 1 as n grows with k fixed, so the first term 
has real part bounded below for n sufficiently large. In particular. 



gin) _ 0{n+l) 



> 1. 



Using the a priori bounds for 6j^ and l/i^ T, we find 



3^"+') and |/ii")|2. 

Hence, 

since kn''^^^^/n~^~^ = 0(721/4-2+0.03+5/3+0.01^ ^ ^^^y ^^ deduce that the second term 
of (jl]) is dominated by kn~^^^~^'^'^ , and then 

^(n+l)| (n)|2 

^(n) _ ^(n+1) ^ ^ 

Changing the value of e appropriately gives the desired result. D 

This lemma is enough for us to estimate the change in ^^" as n grows, and in particular 
to find a limit for the renormalized angle. 

Theorem 7.3. There is a sine-kernel point process {yk)k& such that almost surely, 

^ef =y. + 0((l + A;>-^+^), 

for all n > 1, \k\ < v}l^ and e > 0, where the implied constant may depend on {um)m>i 
and e, but not on n and k. 

Proof. The proof proceeds exactly as in [2J. It is sufficient to prove the result for e equal 
to the inverse of an integer: hence, it is enough to show the estimate for fixed e. By 
symmetry, one can take k > 0. We rearrange the equation in Lemma [7. II to find 

dn) 



(")i^-'^-l)(l + 0(^n-^-)) 



"k 
W|2 _ n('-^-l+2e 



Because almost surely, |/i^ | = 0{n ), we get 

/g(") 

^k 
Using the asymptotic log(l — 5) = —5 + 0(5^) for 5 = o(l), we conclude, if e is small 
enough. 
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Define the random variable L^ = log^^ + X]?=i l/^fc Ti *^ have just shown L^" — 
L]^ = 0{kn^^^^'^) so for k fixed, L^" converges to a hmit L)^ almost surely as n — ?■ oo, 
with |4"^ - 4"^ I = 0(A;n-5+3e)^ ^^^^ 

n— 1 
-in) _ /)W_^ \:^|,,(i)|2 



exp4'^) = 4")exp$:i4^- 



j=k 

(7~i — 1 Tl — 1 



Recall — logn + 'YTj=i ^ = 7 + 0{n ^) where 7 is the Euler-Mascheroni constant. Next 
we define 

71-1 



M«:=E(|^1^'|^-}) 



and observe that each term of the sum is an independent mean-zero random variable. 
Therefore, for k fixed, {M^ )n>k is a martingale. We claim that M^ is bounded in L^; 
in fact. 

It 

so that 

where M];°° is the claimed limit of M^" . To see this, we write 

I (")|2 ei 

where the variables e^ are independent standard exponential random variables. Then we 
compute 

,(«)|2 1^2 _ w / (^- l)ei -62 



E(i/irr--r = E 



n \ n(ei H h e„) 

As shown before in this paper, P(ei + ■ ■ ■ + e„ < ^) = 0{rr'") for all C > 2 so that 

E(l4"^r - ^)' < 0(n-^) + ^E(((n - l)ei - e^ e^f) 



<0(n 



-2\ 



Now, by the triangle inequality and Doob's maximal inequality, for q positive integer, 
k< 2^ 

E(sup(M(~) - Mt^r) < E{{Mt^ - MPr) + E(sup(Mi") - mP)') 

n>2i n>2i 

< EiMf) - Mi^y 

= 0(2-«). 



24 

Hence, 

E 
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r(~) 



sup sup {Ml°°> - Ml 



{n)\2 



< E 



sup sup (Mi'^^-Mi"^ 

A;<2('J+l)/4 29<n<2'J+l 



< E E 

fc<2(9+l)/4 



sup {Mt^-M^-^f 



29<n<2'J+i 



and 



E 



sup sup {m'^^^-MPY 

n>2i k<n^/^ 



r>q 
r>q 



sup sup (M^^ - Ml 

2'-<n<2'-+i A:<ni/4 



3r/4 



0(2 



-3g/4N 



By Markov's inequality, we get 

P(sup sup |Mi°°^ - Mt^\ > 2-^/3) < 22'//3E(sup sup (Mi°°^ - Mi"^) = 0(2-'^/i2) 

n>29 fc<ni/4 n>29 fc<ni/4 



which, by Borel-Cantelh lemma, shows that almost surely for some qo > I, all q > qo, 



n 



> 2" and k < v}l^ satisfy |M^°°^ - Mf\ < 2-«/3. Hence 



|^(oo)_^(n)| 



0{n 



almost surely. Collecting these estimates and applying them to the equation 



'{«) 



3(n) 



r(n)> 



gives us 



exp Vl"' = nOr exp(7 + 0{n-^) + Ml j 



Rearranging, 

n^") = exp(4~) - M^°°) - 7)(1 + 0{kn-'^+^')) =: 27ryfe(l + 0{kn-^^+^'). 

Now, by [2J, {yk)k€Z is a determinantal sine-kernel process, so we have almost surely the 
estimate yk = 0(1 + |A;|), which proves Theorem 17. 3[ D 



8. Weak convergence and renormalization of the eigenvectors 

We are now ready to show that the eigenfunctions /^" of m„ converge in a suitable 
sense. We assume that for a given value of e, the event E from Section [6] holds, which 
happens almost surely. Recall from Theorem 14.21 that we can choose representatives fjf' 
for the eigenvectors of each m„ in such a way that 



(/^r¥/r^^ = E 



/^ 



(n) 



(n) 



where 



fn+l) 



1 J /c 






- f ^ ^ + 



iy„ 



1-A 



(n+l) 



Cn+l 



+ 



kn-l| 



■^ I \(") \("+i)|2 h \(" 



+1)|2' 
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For n = 1, we adopt the convention /} = — ci. We deduce, for n > i, 

n (n) 

3=1 ^j ~ ^k 

The invariance by conjugation of the Haar measures imphes that each eigenvector /^" , 
multiphed by an independent random phase of modulus 1, is a uniform vector on the 
complex sphere S"'~^^. Hence, the scalar product (/^" , e^) converges to zero in proba- 
bility. In order to get a limit which is different from zero, we need to consider a suitable 
normalization. We introduce the following eigenvectors, for n > k: 



(") n(") ^(") 



iJk ■— ^k J I 



k ' 



iW 



where D\. G C is the random variable 



■«-l ^ y(s) _ ^(s + 1) 



or = n(''i'"')*^^^T« 



s=k f^k 

We claim that for each renormalized eigenvector g^\ the scalar product ((7^" , e^) con- 
verges. 

Theorem 8.1. For each k > 1 and i > I, the sequence {{g^ ,e£)}n>kve ^s a martingale 
with respect to the filtration {Bn)n>kve, o,nd the conditional expectation of \{gj^\ei)\^, 
given A, is almost surely hounded when n varies. 

Remark 8.2. More precisely, {{gj^\ei)}n>k\/i is a martingale in the following sense: for 

n > kVi, the conditional expectation of {g^ , e^) given Bn is almost surely well-defined 

and equal to {gj^\ei). However, we do not claim the integrability of {g)^\ee) without 
conditioning. 

Proof. From the equation above, 

n (n) 

{9k .ee}-Dj^ [h^ ) ^ l^-^^ 7(^^\h ' ^^/ 

\{n) x(n+l) n (n) 

^ T^(n) \ - \ Sr^ /^i / An) V 

k („) Z^ ^.{n) .(n+l)^-'j '^^1 

P'k j=l ^j ~ ^k 



(n) _ , (n+1) (n) 

3+k 



Now, recall from Lemma 15^ that the sequence of eigenvectors {fl )n>k is adapted to 
the filtration of cr-algebras {Bn)n>k- If we decompose 

(n) i(n)\ in) I 

then we see that for n fixed, {0j }i<j<n is a family of iid random phases uniformly 

distributed on the unit circle, independent of the cr-algebra Bn- Then since h^ is real 
we have 
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for 1 < A; < n, and we can write 



jt" 



An) 






Wi n(")| fc 



(n) 



A 



(n+l) 



i/^r 



(«)i 



E 







(n) 



l^i 



(«)| 






l<j<n Vk ^V 



(n+l) Jn) 



We would like to compute the conditional expectation of the difference {g'j!''~^'^' — g\, , e^), 
given the a-algebra Bn- We first verify the measurability of each quantity on the right; 
in particular, 



(1) \Dl 



Wl Ai")-Ar^) 






(^7^ — ' ,(n)_ (n+l) are ^measurable. 



(2) <l>^ and (/j , e^) are i3„-measurable. 

(3) {0," (0fc J "'^}i<j<n are iid and independent of Bn- 



JT^k 



We also have the upper bound 

\{9t"-9r,ee)\<\D, 



(n+l) Jn) ^_\| ^ inWl l^fc ~^fc I y^ 



I (")l 



i^r 



(n)| 



An) 



l<j<n I *j 



A "^ - A 



("+i)i 



which shows that ((jf^" — g^^^e^) is uniformly bounded by an ^-measurable quantity 
which is almost surely finite. Therefore we know that the conditional expectation is 
almost surely well-defined. By Fubini's theorem and measurability we get 



E[(^i 



,(") 



9l 



,)\Bn\ = ^V\D 



(")ln(")l^fc 



Ar - A 



(n+l) 



i^r 



(n)| 



E^ 

l<j<n 
3+k 



r 
t' 



Br, 



l/x 



Wi 



in) 



A 



"^J ' inn) V 

(n) ,(n+l)\''i ''^^Z' 



A 



However, by independence. 



E 



(,(") 
r 



Bn 



E 



0} 



0, 



so ((s'l; , e£))n>fcv£ is a (i3„)„>fcvrmartingale. 

To check its conditional boundedness in L^, given A, we need to show 



E[|(^r^e,)n^]+ 5^E[ 



(«+l) 



,(") 



(7r,e,)n^]<oo 



n>fcV£ 



almost surely. The first term smaller than or equal to \\g\. 



(fcV£)||2 



\D 



(fcV£)|2 



which is A- 



measurable and almost surely finite. Hence, it is sufficient to bound the sum. Expanding 

it gives: 



(6) 




E[ 


where 








S = 


■ E 

l<ij<n 
i,j¥'k 



li'^'^-g':\ee)\'\A] = \D 



I \(") _ \("+l)|2 
(n)|2 l^fc ^fc I c 



l/i. 



Wi 



I (")l 



^(n) _ ^(n+1) ^W _ ^{"+1) 



E[ 



t(n) j,(«)/i-(") 



0r(/r,e.)(/r,e,)i^]. 



(n) 
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Now 



E[#Vf \/-^"\e,)(/f ,e,) I A] = E[E[0;"Vf ^ I ^n](i^^"^e,)(/f \e,) | A] 



ein) „ \//(«) 



<5„-E[(/r,e,)(/"\e,)|^] 



where Sij is the Kronecker delta. Thus 



,(")|2 



p) ^= E i,,..)_\,Ui. m<4"'.^«)n-^i- 

l<j<n l-^i ^fc I 

In order to effectively bound this sum we need a posteriori information from the conver- 
gence of the eigenvalues in Section [71 We define the event Fk to be 

F, := {3no > 1, Vn > no, (<\ - 6^"^) A (^i") - 6^"},) > n-^-} 

Lemma 8.3. Fk is A-measurable and holds with probability one. 

Proof. Fk depends only on the eigenangles (hence eigenvalues) and is therefore clearly A- 
measurable. By Theorem 17.31 applied to k — l,k,k + l, we see that each of the associated 
eigenangles satisfies 

ne'f^ =2'Ky^ + 0in--^). 
In particular, 

t^ - C\ = 2^n-\yk - yk-i) + 0{n-i). 
Since {yj}j^z is a sine- kernel point process yk — yk-i ^ 1 almost surely; similar for A; -|- 1 
and k. D 



Now we condition on Fk and estimate the quantities involved in E[| {g)^ Qk i ^^) P 1^]- 
We will begin by estimating D]^ . By Lemma 17.11 almost surely, 

hence 

(8) lA^^) - \t\ = ^r^Vr^P(l + O(n-H^)). 

Similarly, 



|i \("+l)|2 



Now, as in the proof of convergence of the eigenangles, let us consider the set J of indices 



j such that 9^ — vr < 9j < 6'^" -|- vr. From the event F^, we see that for j G J\{k} 
and n large enough, 

|fl(") _ /g("+l)| ^ |/)(") _ zi(n)| _ loi") _ /i("+l)| 



> n-- - o (^c-i/^rra + o(n-t+^))^ 

> n-^-^ - 0(n-2+^) > n-^-\ 
and then 

(9) |Af)-Ar^)|>n-^-^ 

if e is taken small enough. We can get a stronger estimate for \j — k\ sufficiently large. 
In fact, since E^ holds, we have 

(10) |Af^-Ar^)|>|j-A;|n-K 
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These two lower bounds let us estimate 



l/i 



{")|2 



l<j<n l^^j -^k I 



< «,-!+' 



Ei^r-A 



(n) ^(n+l)|-2 

A: I 






< n "'^"'"''(?T,^"'"''n3 + n 3 +''n 3) 



< na 



2 

where we split the sum into the interval with \j — k\ <n^ and its complement. 
Now we can estimate h^ . In fact, 



'^ W|2 



1 + 



1 + 



l^fc ~ ^k I 
\\{s) \(*+l)|2 

Al, At, 



/ 



E 



l/i 



W|2 



\ 



+ 



Iz^., - 1| 



l\W _ \(^+l)|2 h _ \(^+i)|2 

l<j<n l^j '^fc I I-"- ^A: I 



\^^t 



W|2 






d(s+1)n-2. 



l + 0(s-3+^)) 



3(«+l) 



Now, since almost surely, 9^ = 0{l/s), S3+'' = 0{{6f^ ) ^s 3^''), and then 



{^+I)w2,-i+^ 



|\(s) \(«+l)|2 

,{s+i)\j\^_Zj2I^_L - 1 






i/^r 



W|2 



l + 0(s-3 



Now, 



|^{s) _ _)^{^+l)|2 _ iflW fl(s+l)|2 






Using Lemma I7.H one obtains: 



|\(s) \(«+l)|2 



{s+l)\-2|/i(^) /i(s+l)|2 



/ir^^^^^ ^ L = 1 + ^ ^^fc ^fc L ( 1 4. 0(S- 3 



l/^L^ 



(*)|2 



i+i/ii^^r(i+o(.-^- 



,(^)|2 



Applying the bound on |/i^ | given by E2 and changing the value of e gives 



wis) x(^+l)|2 

^(s+1) \^k ~ ^k I 



,(^)|2 



ll^t 



(^)|2 



i+\nr\'+ois-^n 



Thus from the expression 



|^(")|2 _ TT ^(^+i) l-^fc ~\ 



m 



\i^ 



Wi 
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we deduce 

n-l 



\Di'^\' = ll(^ + U'^\' + 0{s-l^n) 



As before, 



s=k 

'n-l -. \ /n-l \ n-l 

exp , J2-]e^vij2U'^\' - M exp5^0(.-|+^) 

.s=k / \s=k / s=k 



n-l 



exp y ^- = k ^nexp[y{l + 0{n ^))] 

s=k 

where 7 is the Euler-Mascheroni constant, 

expE(l4-'l^-i) 

„—l. \ / 



is equal to M„_i — Mk from Section [71 and the last term converges to a limit A^oo with 
error 0(n~3+'^). 
Thus, we have 

|d(")|2 = A;-inexp(7 + M^ - M^ + N^){1 + 0(^-3+^)) 
(11) =:Dkn{l + 0{n--^+')). 

where D^ is a random variable that depends only on k. 

We are now ready to estimate E[|(5f^"' — g^ , e^)^ |^]. In fact we have 

I \(") _ \("+l)|2 

i/"fc 1 

<n\f.twrys 



and 



I (™)|2 

s= E ,,i.)_\L.,, ^ii</!"'-''>i'i-^i 

2 

ieJ,0<|j-/c|<n3 

ieJ,b-fc|>n3 
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Therefore, it is now sufficient to prove that for some e > 0, and almost surely, 

/ 



E 

n>k 



n 



-3+e 



yeJ,o<|j-fc|<n^ 



\ 



+ Y. \3 - kt'n'f^nm 



in) 



A] 



< oo. 



j£ J,|j — fc|>nS 



/ 



It is then sufficient to prove that the expectation of the left-hand side is finite. Since 
J C [k — n — l,k + n], one deduces that it is enough to have 

/ 



E 

n>k 



n 



-3+e 



\o<\j-k\<n'3 



\ 



+ J2 \j - k\-v^'mf^ 



in) 



jgj,n3<|j-fc|<n+l 



< OO. 



Now, since / is, up to a phase of modulus 1, uniform on the sphere S"', one has 



E[\{&\e,)\'] = l/n, 



and then one needs only to check: 



E 



n 



-3+e 



n>k 

which is easy. 



\ 



E ""*' 

\0<\j~k\<n'S 



E b - ^1 



'2^i+^ 



jGJ,n3<|j-fc|<n+l 



< OO, 



D 



Because this martingale is bounded in L^, we have the following immediate corollary. 

Corollary 8.4. Almost surely, for all k G Z and i > 1, the scalar product {g^\ee) 
converges to a limit gt/ when n goes to infinity. 

For each A; G Z, the infinite sequence gk := {gk,i)e>i G C°° can be considered as the 
weak limit of the eigenvector g^ of m„, when n goes to infinity. At the end of the paper, 
we construct a random operator whose eigenvectors are given by the sequences gk, k E Z. 
In order to ensure that these eigenvectors are nontrivial, we need to check that gk is 
almost surely non vanishing. This will be done in the next section. 



9. The law of the coefficients of the eigenvectors of the flow 

In the previous section, we have proven the almost sure convergence of each coordinate 
of the eigenvectors of Un, after normalization by a factor Dj^ . Using the estimate (11 II) . 



we see that almost surely, \D^ \ is equivalent to a nonzero constant times y/n when 
n goes to infinity. It is then more elegant to formulate the result of convergence by 
taking eigenvectors of norm exactly -y/n. Moreover, such a normalization provides the 
distribution of the limiting coordinates of the eigenvectors. The full statement we prove 
is the following: 
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Theorem 9.1. Let {un)n>i be a virtual rotation, following the Haar measure. Fork G Z 
and n > 1, let v^ be a unit eigenvector corresponding to the kth smallest nonnegative 
eigenangle of Un for k > 1, and the (1 — k)th largest strictly negative eigenangle of Un for 
k < 0. Then for all k G Z, there almost surely exist some complex numbers (V'a; )n>i of 
modulus 1, and a sequence {ti:/)e>i, such that for all i > I, 



Almost surely, for all /c e Z, the sequence {tk/)e>i depends, up to a multiplicative factor 
of modulus one, only on the virtual rotation (m„)„>i. Moreover, if {ipk)k£Z is a sequence of 
iid, uniform variables on V, independent of (tk/)e>i, then {il'ktk,i)k&,£>i is an iid family 
of standard complex gaussian variables (EUtpkik/l"^] = 1/ 

Remark 9.2. The vectors f^." are equal to fl , up to a muhiplicative factor of modulus 
1. The independent phases ipk introduced in the last part of the theorem are needed 
in order to get iid complex gaussian variables. This is not the case, for example, if we 
normalize {tk/)e>i in such a way that t^^i G M+. 

Proof. For fixed k E Z, let us first show the existence of the sequence {tk/)i>i. By 
symmetry, one can assume fc > 1: in this case, for all n > k, there exists r^" of modulus 
1 such that fl = Tj^ Vk ■ Hence, 

almost surely converges when n goes to infinity. Now, by the estimate (ITT|1 . \D^ \/^/n 
converges almost surely to a strictly positive constant. One deduces the existence of 
{tk,t)i>i, by taking 

^k ■— ^k ^k ■ 
Let us now check the uniqueness, by supposing that two sequences {tk/)e>i and (t'^ i)e>i 
can be constructed from the same virtual rotation (m„)„>i. In this case, there exist, for 
all n > 1, two unit eigenvectors w^"" and w^ corresponding to the same eigenvalue, and 



such that for all £ > 1, 




(12) 


Vn{w[ \ee) — y tki 

n— >-oo 


and 






71— !>00 



(n) 



Since the eigenvalues are almost surely simple, for all n > 1, there exists x^ ^ U such 

that w)^ = Xk "^k 1 which implies 

(13) xt^^{^t\e,)-^t'. 



k,t 
oo 



By comparing fll2|l and f lT3|) . one deduces that t'^^ = Xktk,e, where Xfc G U denotes the 
limit of any converging subsequence of (xi" )n>i- 

Moreover, let us choose the random vectors {Wk)kez,n>i and the random variables 
{tk,£)k&/>i as measurable functions of (un)n>i, in such a way that f lT2|) is satisfied almost 
surely. Let {4'k)kez be iid random variables, independent of (m„)„>i. For all n > 1, 
the invariance by conjugation of the Haar measure on U{n) implies that the family of 
eigenvectors {ijJkU'k' )i<k<n of Un forms a Haar-distributed unitary matrix in U{n). One 
deduces that if L is a finite set of strictly positive integers, and if i^ is a finite set of 
integers, then 
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converges in law to a family of iid standard complex gaussian variables. Since this fam- 
ily of variables also converges almost surely, one deduces that the limiting variables 
{ipk'tk/)ki^K/&L are iid standard complex and gaussian. Since the finite sets K and L can 
be taken arbitrarily, we are done. D 

One deduces immediately the following: 

Corollary 9.3. The limiting coordinates Qk/ introduced in Corollary \8.4\ are almost surely 
different from zero. 

Proof. We know that for a suitable normalization of t^/, 

^.,= lim|Di")|($i"Vr,e,), 

tk,e= limv/^($i"Vf\e,), 
and by ([II]), 

^/D~k= lim|Di"Vv^- 

71— >00 

Combining these limits gives 

(14) gk,e = \/D'ktk/, 

which is almost surely nonzero, since tk^i is a standard complex gaussian variable. D 

The eigenspaces of «„, generated by the vectors /^" , can also be considered as ele- 
ments of the projective space P'^~^(C). Moreover, one can define the infinite-dimensional 
projective space P°°(C), as the space of nonzero infinite sequences of complex numbers, 
quotiented by scalar multiplication. The convergence of renormalized eigenvectors proven 
above can be viewed as a convergence of the corresponding points on the projective spaces. 

There exists a uniform measure on all these projective spaces, obtained by taking the 
equivalence class of a sequence of iid standard complex gaussian variables. For n > 1 
finite, the uniform measure on P'^(C) can also be obtained from a uniform point on the 
sphere in C"^"*^. For m < n E NU {oo}, there exists a natural projection Il^.m from P"(C) 
to P™(C), obtained by taking only the m + 1 first coordinates of the sequences, and 
this projection is well-defined when these coordinates are not all vanishing: in particular, 
almost surely under the uniform measure on P"(C). Note that the image of this measure 
by n„,m is the uniform measure on P™'(C). Moreover, we can define the notion of weak 
convergence on the projective spaces as follows. Let (x„)„>i be a sequence such that 
Xn e P"(C) for all n > 1 and let x^o ^ P°°(C). We say that (x„)„>i weakly converges to 
Xoo if and only if the following holds: for all m > 1 such that the m + 1 first coordinates 
of Xoo are not all vanishing, the projection Yin,m,{xn) € P'^(C) is well-defined for n large 
enough and tends to Iloo,m{,Xoo) when n goes to infinity. From the previous result, we can 
easily deduce the following: 

Theorem 9.4. Let («„)„>! he a virtual rotation, following the Haar measure. For k eT, 
and n > 1, let x^ E P"^-'^(C) be the eigenspace corresponding to the kth smallest non- 
negative eigenangle of u^ for k >1, and the (1 — k)th largest strictly negative eigenangle 
of Un for k < (this eigenspace is almost surely one- dimensional). Then, there almost 
surely exists some random points {x^ )kez in P°°(C) such that for all k E Z, x^ weakly 
converges to x)^ when n goes to infinity. The points {x^ )k<^'L are represented by the 
sequences {tk/)e>i, k E Tj given above: they are independent and uniform on P°°(C). 
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10. A FLOW OF OPERATORS ON A RANDOM SPACE 

For each a G M, let (q;„)„>i be a sequence such that a„ is equivalent to an when n 
goes to infinity. For n > 1, A; G Z, we have 

„,ar,An) _ je['''>a„ An) 
"'n Jk ~ ^ Jk ■ 

Now, e*^fe "" tends to e^*'^"^* and after normalization, the coordinates of f]^ tend to the 
corresponding coordinates of the sequence {tk/)i,>i. It is then natural to expect that, in 
a sense which needs to be made precise, u"" tends to some operator f/, acting on some 
infinite sequences, such that 

This motivates the following definition: 

Definition 10.1. The space £ is the random vector subspace of C°°, (algebraically) 

generated by the sequences {tk/)e.>i, or equivalently, {gk,e)e>i, for A; G Z. 

For a G M, the operator f/" is the unique linear application from £ to S such that for all 

k ez, 

or equivalently. 

Remark 10.2. For each k, the sequence {tk/)e>i is almost surely well-defined up to a 
multiplicative constant of modulus 1. Moreover, the sequences {tk/)e>i for fc G Z are 
a.s. linearly independent, since for a suitable normalization, {tk/)k&z,e>i are iid standard 
complex gaussian. This ensures that the definition of f/" given above is meaningful. The 
notation [/" is motivated by the immediate fact that (f/")„giR is a fiow of operators on 
S, i.e. f/° = Is and f/"+'' = f/^f/^^ for all a,/3eR. 

As suggested before, we expect that f/" is a kind of limit for u"" when n goes to infinity. 
Of course, these operators do not act on the same space, so we need to be more precise. 

Theorem 10.3. Almost surely, for any sequence {si)e>i in £ and for all integers m>l, 

K"((s.)l<.<n)]^-^[f/"((s.).>l)]^, 
n— )-oo 

where [■]m denotes the mth coordinate of a vector or a sequence. 

By linearity, it is sufficient to show the theorem for (s£)^>i = {gk/)(.>i- Hence, Theorem 
110.31 can be deduced from the following proposition: 

Proposition 10.4. For all k E 1^, i > 1, one has almost surely: 

as n^ oo, for gk[n] := {gk/)i<£<n- 

In the next section, we will give a more intrinsic way to define a fiow of operators 
similar to {U"')a£R- In order to make this construction, we need a more precise and 
stronger result than Proposition 110. 4[ which is given by the two following propositions: 

Proposition 10.5. Let e > 0. Almost surely, for all k E Z, we have the following. 

• \\gk['n]\\ is equivalent to a strictly positive random variable times -y/n, when n goes 
to infinity. 

• \\9k[n] -S'i"^!! = 0,(n^+'). 
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• For any T > and 6 G (0, 1/6), 

sup sup \\uT9k[n] - e^^^'^^fcMII = 0{n^-^). 

ae[-T,T] a,jG[n(a-n-*),n(«+n-'')] 

Proposition 10.6. Almost surely, for all k & Z, i > 1, a,'y & M, and for all sequences 
(«n)n>i (ind (7n)n>i such that a„/n = a + o{n~^) and jn/n = 7 + o{n~^) for some 
6 e [0,1/6), 

{K"{g,[n]) - e'^^'^y''g,[n],ul"{e,)) = o{n-'), 
when n goes to infinity. Moreover, for 6 G (0, 1/6), we get the uniform estimate: 

sup K"igk[n]) - e2™^^(7,[n],<"(e,)) = 0(n-^). 

oi„£[n{a—n ^),n(a+n *)] 
7ne[n(7-n-''),n(7+n-^)] 

Proof. We now prove Propositions 110.51 and 110.61 it is clear that this last proposition 
implies Proposition 110.41 (by taking 6 = •jn = 0) . 
Using (IT^ . we get the following: 

n 

ll^ffcNlP = \Dk\'^\tk,e\'^ 
e=i 

where {tk/)e>i are iid standard complex gaussian variables. By the law of large numbers, 
[[(^^[n]!^ is a.s. equivalent to n\Dk\ when n goes to infinity, which shows the first item of 
Proposition 110.51 

Now, the third item can quickly be deduced from the second one. Indeed, if we have 
the estimate H^ffcH — dk W ~ 0(n3+^) for all e > 0, then, since 5 < 1/6, we also have 
llfi'fcW "fl'fc II = 0(^2^'^), which implies, for all a G [—T,T] and a„ G [n{a — n~^),n{a + 
n-% 

IK-^fcM-e^^^'^^fcMII 



< ll<"(^;.M - ^f )ll + IK-^f - e^™^''^^ II + ||e^™^H^4^] - ^i"^^ 
= 2||^,M-^f|| + |e-"^i"'-e^™^MII^i"^ll- 



Now, 



and 



lkN-^fll = oM-^), 
l^f II < II^.NII + II^.M - ^i")|| = o(v^) + oM"^) = oiV^), 






< \9^f^^\\an - an\ + \a\\n9l^''^ - 2nyk\ 

= 0{l/n)0{n^-^) + 0{n-'^*) = 0{n-^), 

where the implied constant does not depend on a and «„ (recall that |a| is assumed to 
be uniformly bounded by T). Note that we used Theorem 17.31 (for k fixed and e = 1/12) 
at the last step of the computation. This gives: 

\\Krg,[n]-e'-'-y^gk[n]\\=0in'2-'), 

uniformly with respect to a and «„, i.e. the third item of Proposition 110.51 In order to 
complete the proof of this proposition, it then remains to show the second item, which 
needs several intermediate steps. The proof of this convergence needs several steps. 
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Lemma 10.7. For fixed k > 1, e > 0, there exists a A-measurahle random variable 
M > 0, such that, almost surely, for all N > n > k , 



Sn) 



nWar - 9k[n]\\' \A] < M 



na 



and 



» J^)\\2 



m9r-9n\'\A]<M{N 



n)n 3 



h+e 



^N) 



where g),^ denotes the vector obtained by taking the n first coordinates of g)^ 
Proof. One has, for n > 1, and A^ G {n, ra + 1, ra + 2, ...} U {c)o}, 

n 



(TV) 



S^) 



(m) 



f*^^ 9kn '■~ fl'fcN- By the martingale property satisfied by the scalar products {g^ ,ei) 
for m > k, 

m9t'\e,)-{gf\e,)\'\A]= J^ Af^ 

n<m<N 



where 



A^r^=E[\{gt^'\e,)-{g^-\ee)\'\Al 



and then by IQ and ([7]), 






\{m) \{m+l)|2 



(m) I 



E 



\fi 



{m)|2 



(m) , (m+1) I 



-E 



\h^k I l<j<m,j^k l^j '^fc I 



p(m) 



l(/r'e^)rl-4 



(15) 



<Mim-3+^ ^ |aJ™)-A 



(m) -i (m+1) I -2 



E 



l(/^\e,)n^" 



ieJm\{fc} 



where J^ is the random set of m consecutive integers, such that 9\. 
^k + ^) ^iid where 



(-+1) _ vr < 0(™) < 



1 := sup m \Dl '^ 



Am) x(m+l)|2| (m)|2 



''T^^k,l<j<m,jj^k 



I ("*)|9 



From the estimates (|H]), f lTT]) . the fact that \fi- \ and |yU^ | are almost surely dominated 



(m)|2 



,M|2 



by m ^+2 and the bound 6'^™ = 0(l/m), one deduces that the ^-measurable quantity 
Ml is almost surely finite. Similarly, one has, for all j G Jm\{k}, 



(16) 
where 



|^M_^(„^+l)|-2<^ / 2+e/^ 



^ 3 -t-t 



<Mo 



ms" 



I^-J'l 



Mo:-- 



sup |A^- — A 

m>fc,jeJm\{A;} 



(m) A (m+1) I -2 I ,^2+ 



m""^' A 



10 

m 3 ' 



l^-il 



Now, M2 is ^-measurable and, by the estimates (jH]) and flTU]) . it is almost surely finite. 
From flT^ and fITB]) . we deduce (with the change of variable p = j — k): 



aJ"^)< MiM2m-5+2^ y ^E 

\p\ 

pG{-m-l,-m+l,...,-l,l,...m+l} ' ' 



.(m) 



\{nZee)['\A 
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Therefore, 

1 



n{9r.ed-{9T\e,)\'\A] 

< M1M2 Y^ 



m-H2^ 



^ \p\ 



n<m<N pe{— m,— m+1,..., — l,l,...m} 

and then, summing for C between 1 and n. 



\{ftle,)?\A 



< 



M1M2 Yl 



m-5+2^ 



E 



n<m<N pG{—m,— m+1,..., — 1,1,... m} 

Let us first suppose that A^ is finite. One has 

n m 

which imphes. 



bl 



E 



Ek/S^^)|'i-4 



c{»")||2 



\\Jk+p\\ — -■-' 



m 3+2^ 



E 

pS{— m,— m+l,..., — l,l,...7?i} 

< M3 ^ m-^+3^ < Ms{N - m)m-^^^' 



n<m<N 



1 

H 



where 



n<m<N 



M3 := M1M2 sup I m"" V — 

pG|— m,—m+l,...,~ 1,1,... m| 



< 00 



is ^-measurable. Hence, we have the desired bound, after changing e and taking M = M3. 
In the case where N is infinite, we use Lemma [6.31 to get the estimate 



E 



P(m) 



\{nz^^)n^ 



< MiTJl' 



-1+e 



where 



M4 := sup m^~'E 

m>l,l<k,i<m 

is v4-measurable and almost surely finite. Hence, 

< M1M2M4 Y ^~^' 



\{ft\ee)\'\A 



+3e 



E 



n 



m>n. pG{— m.,-m+l,..., — l,l,...m} 

which is easily dominated by a finite, ^-measurable quantity, multiplied by 713+^*^. D 

We have now an L^ bound on \\g^ — QknH, conditionally on A. The next goal is 
to deduce an almost sure bound, by using Borel-Cantelli lemma. This cannot be made 
directly, since the corresponding probabilities do not decay sufficiently fast, but one can 
solve this problem by using subsequences. 

Lemma 10.8. For fixed k > 1, e > 0, let u := 1 + [3/eJ, and for r > 1, let Ur := k + r'^ . 
Then, almost surely, 
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Proof. From Leinina riO. 71 (applied to e/2 instead of e), there exists M' almost surely finite 
and ^-measurable such that for all r > 1, 



and then 



M' + l 



E 






< 



M' 



M' + 



ivi ^-^ 



-'/^ < 3. 



r>l 



By taking the expectation and using the fact that 1/(M'+1) is ^-measurable, one deduces 



E 



— El 

' + 1 ^ 

r>l 



M' + 



\\gl;'''^-gklnrW>n?' 



<3, 



and then 



f' + i ^ iifli"'-' 



M' + 1 ■^ \\gr'-9k[nrW>n. 



almost surely. Hence, almost surely, \\gl. 
r > 1. 



(nr) 



gk[nr 



2 ,, < OO 



< nf- 



for all but finitely many 

D 



The next lemma gives a way to go from a given subsequence to a less sparse subse- 
quence. We will say that a value 5 G (0, 1] is good if the conclusion of Lemma nO.SI remains 
valid after replacing the sequence (k + r'^)r>i by {k + [r^^^\)r>i, the brackets denoting 
the integer part: from Lemma ri0.8[ we know that 1/z/ is good. 

Lemma 10.9. If S & \^l'^i iy ~ l)/'^] ^■5 good, then b + 1/z/ is good. 

Proof. For r > 1, let Ur := k+ \^r^/(^+^/'^)^ and let A^^ the smallest element of the sequence 
{k + [s^^^\)s>i, such that A^^ > ""r- The sequence {Nr)r>i, as {nr)r>i, tends to infinity 
with r: moreover, it is a subsequence of {k + [s^^^\)s>i (but some terms of this sequence 
may appear several times in {Nr)r>i)- By assumption, one has almost surely: 

\\gi'''^-g,[NrW = 0{N}^'). 

Now, it is easy to check that N,. = 0{n^) (with a constant depending only on 5). Hence, 
by restricting the vectors to their n^ first coordinates, one deduces, almost surely. 



Oin'r 



(17) \\9t!;y9k[nrW<\\9f-^-gk[Nr]\\ 

On the other hand, the distance between two consecutive terms of the sequence {k + 
[s"'^/'^J)s>i satisfies the following: 

{k + [{s + lf"\ )-{k+ [s^"\ ) = {s + lf"- s^" + 0(1) 
which implies 



Nr — rij. = 0(n. 



1-S\ 



where the implied constant depends only on k and 6. One deduces, by Lemma [10.71 

2 . ^ _5 



(nr) 



,(A'r)||2 



n\\9r - 9r;x \A] < M' {Nr - Ur 



[rir 



'*+5 <M'nl^^' 
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and then 

where the impled constant depends only on k and 6. Since the exponent of r is strictly 
smaller than —1, one deduces, by using Borel-Cantelli lemma similarly as in the proof of 
Lemma [10.81 that almost surely, 

(18) \\9':^-^-a\\' = 0inl^'). 

Combining (TT7|) and (TT8l) gives the desired result. 

D 

By applying Lemma 110.81 and {u — 1 times) Lemma 110. 9[ one deduces that 1 is good, 
which gives the second item of Proposition 110.51 for k > 1: the situation for fc < is 
similar. 

Let us now prove Proposition 110.61 By triangle inequality, we get: 

\K"{g,[n])-e'-^-y>=g,[n],ul-{ee))\ 

< \{K"i9k[n] - ^i"^),«I"(e,))| + |«"^f - e^™^'=^i"\«;^"(e,))| 

+ Ke^^^H^i"^ -^4^]),<"(e,))| < K^f -^.M,<"-""(e.))| 

+ |e^""'^"' -e^™^M|(^i"\tx;^"(e,))| + {{gt"^ - 9kHul"{e,))\. 

To prove the first part or the proposition, since the sequence (a„ — 7„)„>i satisfies exactly 
the same assumptions as (7n)„>i (replacing 7 by 7 — a), it is sufficient to prove the almost 
sure estimates: 

(19) |e^°»^i"' - e^™^'=|l(^i"U;i"(e£))| = o(n-^) 
and 

(20) \{glr^-g,[nUl-{ee))\=o{n~'). 

For the uniform part of the Proposition 110.61 it is sufficient to prove the same estimates, 
with o{n~^) replaced by 0{n~^), this bound being uniform with respect to a„ G [n{a — 
n~^), n{a + n~^)], and 7„ G [^(7 — 2n~'^), 71(7 + 2n~^)] (the factor 2 comes from the term 
where 7^ plays the role of 7^ — a„). Now, under the assumption of the first part of the 
proposition, fITI?]) is a consequence of the two following estimates: 

= 0{l/n)o{n^-^) + 0{n''^^) = o{n-^), 



and 



.(") ,,7nY^ \\| _ |/„-7n^(") , \| _ |/^-n.4"' W 



= \{gt\ee)\=0{l), 

the last estimate coming from the almost sure convergence of {g^ , e^) towards gk/- This 
computation is also available (after replacing small o by big O) for the second part of 
the proposition, since all the estimates are easily checked to be uniform with respect to 
an G [n{a — n~^),n{a + n~^)] and 7„ G [77.(7 ~ 2n~^),n{'~f + 2n~^)]. 

It remains to prove fl20l) . From the second item of Proposition 110.51 and the fact that 
6 < 1/6 we have the estimate: 

\\9t'^'9k[n]\\=o{n^-'') 
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for some 5' > S. Hence, in order to complete the proof of Proposition I1U.6[ is sufficient 
to show the following property of delocalization: 

(21) sup \i9':''-9^[nUnee))\ ^ o{n-'^^^'-% 

7ne[n{7-2n-«),n(7+2n-«)] \\g]^' - gk[n\\\ 

This will be a consequence of the following lemma: 
Lemma 10.10. For all n > 1, the quotient 



is a P random variable of parameters 1 andn — 1, i.e. it has the same law as \xi\^ , where 
(xi, . . . , Xn) is a uniform vector on the complex sphere S^ . 

Proof. Let m > 1, let a be a random matrix in U{m), independent of («„)«>! and 
following the Haar measure. Let us define (tij^)„>i as the unique virtual isometry such 
that for all n > m. 



, _i o \ a 



The invariance by conjugation of the Haar measure on the space of virtual isometrics 
implies that (m^)„>i has the same law as {un)n>i- Let {g), )„>i be the sequence of 
eigenvectors constructed from (u^)„>i in the same way as {g)^ )n>i is constructed from 
{un)n>i- Since (M^)n>i has the same law as («„)„>!, one deduces that almost surely, each 
coordinate of g^ with a given index converges to a limit when n goes to infinity: let g'j^ 
be the corresponding limiting sequence, which is the analog of gu when {un)n>i is replaced 
by {u'^)n>i- One knows that for n > m, 

si"' - ( / ° ) si"' 

y "J ^n—m J 

is an eigenvector of (m^)„>i, corresponding to the eigenvalue A^" (recall that for n > m, 
Un and u'^ have the same eigenvalues). Hence, almost surely, there exists Hn ^ C* such 
that g)^ = Undk ■ Moreover, from ([5]), and from the fact that gj!^ is orthogonal to /j" 
for j 7^ k, one obtains that 

and then 

; ~(n+l) ~{n) ~(")\ n 

Kl^n+lQk^n ~ '^nQk i '^nQk I — U- 

Applying diag(cr^^, In-m) to the vectors in the scalar product and dividing by |k„P gives: 

(n+l) \ 

^^+1 3k,n in) in) \ _ „ 

iJk ^iJk ~ ^^ 



Kn 



and then 




C -^r,^n=o, 



which implies fi;„+i = 

(9t:'\9t') = {9t:''-9':\9':') + ii^rir = ii^rir > o. 
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Hence, one has (7^" = Hmgf!^ for all n > m: by taking the limit for n — )■ 00, one deduces 
that g[ = t^rnQk, where gk is the infinite sequence obtained from gk by applying a to the 
vector formed by the m first coordinates, and by letting the coordinates fixed for the 
indices strictly larger than m. One deduces the following (with obvious notation): 

(^f^' - 9'Aml {uLr-{ee,)) = K^ {-g'r^ - h[m], K)^-(e,J) 

= i^m {crig'r^) - cr{gk[m]),a{u^y"'a~^{eeo)) 
= Km {gt^^ - gk[m], (Mm)^'"o-"^(ef J) 
= Km {{umy^"'{g'r^ - gklm]), (T~^(e<.J). 
Similarly, 

II (™)' n 1I|2 I |2||~(m.) ~r ]||2 I |2|| / {m)\ / r im|2 

\\9k -^feNII =\Km\ Wgl - gk[m\\\ =\Km\ \\(T{gl ') - a{gk[m\)\\ 

= I'^mHls'i™^ - gk[m]\\^ = \Km\'^\\{umy"'{gk'^^ -gk[m])\\^. 



Hence 



(22) 



\{9t^'-9'kH,iu'mr-ieJ)\' \{x,y)\' 



119^^^' - 9'k[mW INI' 



where x = {um)^'^"^{gk ~ 9k\f^) and y = a~^{e(,f^) is independent of x (since a is 
independent of {Un)n>i) and uniform on the unit sphere of C"* (since a is uniform on 
U{m)). One deduces that the left-hand side of (12 2 p is a beta random variable with 
parameters 1 and m — 1. Since {u'^)n>i and («„)„>! have the same distribution, one can 
remove the primes from this left-hand side, which gives the announced result. D 

Now, applying Lemma flU-lOt we get: 

P ( sup \{9'^''9k[nWnie,))\ ^ ^^...S 

\7ne[n{7-2n-*),n{7+2n-*)] | \g]^' - gk [u] \ \ 

< (1 + An~^)F[P{l, n-l)> n-'+^^^'~^^] = O f (1 + An,-^) e-"'^''""/6 



which by Borel-Cantelli lemma, gives the estimate fl2T|) . and then completes the proof of 
Proposition 110.61 D 

11. A MORE INTRINSIC DEFINITION OF THE LIMITING OPERATOR 

The random space S and the fiow (t/")aeM of random operators on S defined previously 
have the disadvantage to involve explicitly the limiting eigenvectors of {un)n>i- This 
definition is artificial. In this section, we construct a random space of sequences which 
contains S and a fiow of operators which restricts to U"' on £. This random space is not 
constructed directly from the eigenvectors, but is rather the space of sequences on which 
the action of the finite matrices (wn).^! converges in a suitable way. 

All probabilistic statements in this section are assumed to hold almost surely. For this 
section, let B{x,r) C M denote the closed ball in R with radius r. 

We define our random space as follows. 

Definition 11.1. Let J-" denote the space of sequences {we)e>i such that for all a G M, 
there exists a sequence V^w, satisfying the following properties {v[n] denoting, for any 
sequence v, the vector {vi)i<e<n G C" for all n > 1): 
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• For all f > 1, a, 7 G M and < 6' < 6 < 1/6, there is a constant C > such that 

_sup |«"(«;[n]) - {V'^w)[n],ul"{ee))\ < Cn-^' . 

an£B(an,n?-^ ) 
7„eS{7n,nl-*) 

• For all T > and < 5' < 5 < 1/6, there is a constant C > such that 

sup sup IK-wN - (V^"w)[n]|| < C (n^'^'\ . 

Here the constants may depend on the choices of £, a, 7, 5', 5, and T. 

Note that the first condition implies, taking a^ = \an\ and 7„ = 0, that for w E J-', 
{V°'w)£ is the limit of {un {w[n]) , ei) when n goes to infinity. Hence, V"w is uniquely 
determined. 

The above definition gives the random space J-' and the family of operators V"' together. 
It turns out that this definition gives a flow of operators on a vector space which restricts 
to S in the natural way. 

Theorem 11.2. The set J-' is a vector space. There is a family {V°')ainR of linear maps, 
V"' : J-" — 7> J-", given by the correspondence w h- )■ V^w. The family satisfies the semigroup 
properties V'^ = id and V^-V^ = V^^^ for all a, /3 G M. Moreover, almost surely for all 
A; G Z and a E^, one has {gk/)e>i G J-" and V°'{{gk/)e>i) = {e'^'^^°'^''gk/)e>i, so that S is 
a subspace of T and that f/" is the restriction ofV" to S. 

Proof. We begin by showing that J-" is a vector space. Clearly it suffices, for Wi,W2 E J^ 
and A G C, to show that \wi + W2 E J^. Let w = \wi + W2 as sequences and define 
Wa = XV^Wi + V"w2. Now consider all £ > 1, a;,7 G M and < 6 < 1/6. For any 
an G B{an,n^~^) and 7„ G B{jn,n^~^), by the triangle inequality 

\{ul-{w[n])-wM.ul-{e,))\ < |A||«"KM) - {V^w,)[n],ul-{e,))\ 

+ \{K-{w2[n]) - {V^W2)[niul-{e,))l 

and 

\\ui-w[n]-wM\\ < \M ll<"^iN - (v^"wi)MII + ll<"^^2M - (V"W2)M||. 

We deduce that w E J^ with V^w = Wa and so J-" is a vector space as required. 

To see that V'^ = id, it suffices to note that {V^w)e is the limit of {un {w[n]),ei) = 
{w, Ci) as n goes to 00. 

Now, let w e J^, a,(3 eR, and let us show that V^w G J" and V^iV^w) = V^+'^w, 
which proves the stability of J-' by V^ and the group property. The linearity of V, of 
course, follows by construction. We need to show: 

• For all£ > 1, a, 7 G M and < 5 < 5' < 1/6, 

_sup \{u'^"{V^w)[n]) - {V''^^w)[n],ul-{ee))\ = Oin'^'). 

an&B{an,n^^ ) 
7„gB{7n,ni-*) 

• For all T > and < 5 < 5' < 1/6, 

sup sup I \u^" (V'^w) [n] - (y"+^w) [n]\\=0 (n^^-^') . 
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Let US therefore fix the parameters i > 1, T > 0, a E [— T, T], 7 G M, and < S' < S < 
1/6. Consider any choice of «„ G B{an,n^~^) and 7„ G B{'yn,n^~^). Define the sequence 

(/3n)n>i by 

_ f [/3nJ , a„ > an 

It is easy to check that /3„ G B{f3n, n^~^) and «„ + /3„ G -B((q; + /3)n, n^~^). Now, by the 
triangle inequahty, 

\{u:-iV^w)[n] - iV^+^w)[n],ul"iee))\ < |«"(\/^^)N - <"+^"(ti;[n]),<"(e,))| 

+ |«"+^"(w[n]) - {V''+^w)[n],ul"{ee))\. 
By unitary invariance, we see that 

|«"(\/^^)M - (\/"+^^)[n],<"(e,))| < |((\/'^t/;)M - w^(«;M),w;^"-""(e,))| 

+ |«"+'5"(u^M) - (\/"+^u;)M,<"(e,))| 

and therefore, taking supremums, 

_sup \{u'^"iV^w)[n] - iV''^^w)[n],ul"iee))\ 

7nG"B(7n,ni-'^) 

< _sup |((\/^w)[n]-M^(w[n]),<"-""(e^))| 

aneB{an,n^'^) 
7neB(7n,nl-*) 

+ _sup I (m^"+^" (w [n] ) - (K"+^w;) [n] , m;^" (e,) ) | . 

a„G-B(an,rai"*) 
7nei?(7n,nl-'') 

Denoting a^ = «„ + /3„ and 7^ = 7n — ctn, we have |/3„ — n/3| < 1, |q;^ — n(Q; + /3)\ < n^~^ 
and |7^ — n(7 — q;)| < 2n^~'^. Hence, 



_sup \{u'^-{V^w)[n] - (V^"+'^w)[n],<"(e^))| 

7neB(7n,ni-'^) 

< s_up \{{V^w)[n]-ut{w[n]),u^J^{ee))\ 

_P„eB{[Sn,l) 
7;GS((7-a)n,2nl-*) 

+ _ sup \{u<{w[n])-{V^^%)Hul-{ee))\ 

a'„eB{{a+fS)n,n'^-^) 
7ne"B(7ra,ni-'') 



S ^ _.")^ 



For 6" G ((5', (5), we deduce, for n large enough (so that 2n ^ < n 
_sup |«"(\/'^«;)[n] - {V''+f^w)[n],ul"{ee))\ 

an&B(an,n^^^) 
7„eB(7n,ni-'^) 

< _sup \{iV^w)[n]-u^^"iw[n]),u'^-iee))\ 

l3„&B{l3n,n^-^") 
7;GS({7-a)n,ni-^") 

+ _ sup \{u<{w[n])-{V"+^w)[n],ul-{ee))\ 

a'„eB{{a+fS)n,n'^-^) 
■yneB{-yn,n^-^) 
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Now, this last expression is dominated by n'^ , by definition of tfie space J-' and tlie maps 
V^ and V°'~^^. By a similar computation for the second part of the definition of J-', we 
get 

\\u'^"{V^w)[n] - {V''+^w)[n]\\ < \\u^"{V^w)[n] - u'^-+^"w[n]\\ 

+ ||<"+^"w[n]-(r"+^w;)[n]|| 

which, again by unitary invariance, reduces to 

\K-{V'^w)[n] - (r"+^w)[n]|| < \\{V^w)[n] - u^^-w[n]\\ + |K"+^"wH - (\/°+^w)[n]|| 

Defining a' := a + (3, a'^ := q;„ + /3„, T' := T + \(3\ and taking the supremum for 
a G [—T,T] and a„ G B{an,n^~^), we get 

sup _sup I \u^" {V%) [n] - (V+'^w) [n] \ \ 

< sup \\iV'^w)[n]-u'^"w[n]\\ 

+ sup sup ||M""w[n] — (V" w)[n]|| 

Both of the quantities on the right hand side are bounded by 0{n2~^ ). We have therefore 
proven the stability of J-" by the family of maps {V")a(^M. and the semigroup property. 

It remains to show that the eigenvectors are contained in J-" and V^ restricts to f/" 
there, but this is an immediate consequence of Propositions I1U.5I and I1U.6I D 

We have now defined a random space J-" containing S, on which the flow {U")a(=:R is 
extended to a flow of operators (y"')aeR- The definition of J-" is quite complicated, but it 
has the strong advantage, compared with the case of £^, to be given intrinsically in terms 
of {un)n>i, without referring explicitly to eigenvectors. A natural question which can now 
be asked is the following: by extending the space £ to J-', do there exist eigenvectors of 
the flow which are not contained in £7 The answer is negative, in the following precise 
sense. 

Definition 11.3. For w E T different from zero, we say that w is an eigenvector of the 
flow {V°')aeR, if and only if there exists x ^ ^^ such that ^"(u;) = e'^'^'^^w for all a; G M. 

Using the following definition, we see that for all A; G Z, the sequences (tfc,^)^>i are 
eigenvectors of (V")agR, corresponding to x = Vk- The following result shows that these 
sequences are essentially the only eigenvectors of the flow: 

Theorem 11.4. The only eigenvectors of {V°')aeK O'l"^ the non-zero sequences which are 
proportional to {tk/)e>i (or {gk,e)e>i) for some A; G Z. 

From this theorem, we deduce in particular that the set of parameters x associated to 
the flow {V'^)a£R is a determinantal sine-kernel process. 

Proof. One important ingredient in the proof is the following lemma, which shows that 
a non-zero element in J-" cannot be too small, in a sense which can be made precise. 

Lemma 11.5. Let w be an element of J^ , such that there exists 5 > 0, and a strictly 

increasing sequence {nq)q>i of integers, for which \\w[nq]\\ = Oijiq ). Then, w is equal 
to zero. 
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Proof. In the first part of tlie definition of J-", let us take (5/3) A (1/7) instead of S, and 
a = 7„ = 0. Since V^w = w, we deduce, for all i > I, 

sup \{u'^"{w[n])-w[n],ei)\ = 0{n-^'), 

\a„\<n^-^ 

for (say) S' = (6/4) A (1/8) > 0. Hence, if we decompose w[n] in terms of eigenvectors of 



itrt.' 



w 



[n]=j:€'ft\ 



k=l 



we get: 



sup 



Wi 



J2(^t'r"vt'\ft\ee) 



fc=i 



0(n- 



Taking p = [n^ sj + 1 and averaging for q;„ G {0, 1, . . . ,p — 1} gives: 



where 



n 

n— >oo ' ' 
k=l 






Moreover, for all n > 1 and k G {l,...,n}, K/^" , e^)p is a beta random variable of 
parameters 1 and n — 1. From the Borel-Cantelli lemma, we deduce that almost surely. 



(/f,e,) = 0(n-M), 



and then, using Cauchy-Schwarz inequality, we get: 



,A;=1 



\we\ ^ n 2^3 



1/2 / „ \ 1/2 

(")m2 \ ( Y^ |ri(")|2 

vA:=l 

1/2 



(n)M2 



(23) =n-^+3||w;[n]||K^|Mp(A 

\fc=i 

Let us now assume the following bound 

n 

(24) El^^(^l"^)l' = 0(^') 



fc=i 



Then, for ^ = n^ (g > 1), the right-hand side of (123!) is dominated by n^ , which is 

only possible if W£ = 0. Hence, Lemma [11.51 is proven if we are able to show (^^. Now, 

we have 

p-i 



which implies 



j=-p+i 



n ^ p-1 

7 1 ^ ^ — ^ I 1 



fc=l 



i=~p+i 
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By [S], it is known, that for all integers ji,J2, ■ ■ ■ ,jr such that |ji| + IJ2I + ■ ■ ■ + \jr\ < n, 
one has 



E 



riTrl 



ut 



E 



n^i 



s=l 



where Zq = n, (^j/a/J)j>i are iid standard complex gaussian variables (E[|Zjp] = j), and 
Z_j = Zj for j > 1. For a fixed integer A > 0, we deduce, from the fact that p = o{n), 
that for n large enough: 

s A' 



E 



^(n)-M2 



.fc=l 



E[W' 



where 



^-^ E ip-\mr 



j=-p+i 

The variable ly is a real-valued gaussian variable, with expectation n/p (coming from 
the term j = 0), and variance 

^ E wp-\Ar<^ E p'^i- 

i<|j|<p-i i<lil<P-i 

Hence 

IIH^IIl^ < (n/p) + \\W-{n/p)\\LA < {n/p) + ||iV(0, 1)|U.4 < c{A)+n/p, 

where c{A) > depends only on A. We deduce (for A fixed): 

xA' 



E 



,W„2 



vfc=l 



< {n/p)^ < n^'l\ 



Using Markov's inequality, we get: 



P 



Y.\M,{\t^)\^>n' 



k=\ 



< n-^^^ 



E \^v^< 



(")\|2 



.fc=l 



< „-2A5/3 



Taking any A strictly larger that 3/25 and using the Borel-Cantelli lemma gives the 
estimate (El). □ 



Now, let w be an eigenvector of (V^")oeK, and let x be the corresponding eigenvalue. 
Taking, in the second part of the definition of J-", T = 1, a = j/n for j G {0, 1, . . . , n — 1}, 
a„ = an = j, we obtain, for any b G (0, 1/6): 



or equivalently, 
where 



sup \\uiw\n\ - e'^'^^^'/^wWII =0[nh-^ 

0<i<n-l 



w\n\ = <e-2^"»/"w[n] + w„,j. 



sup 

0<j<n- 



\\vnj\\ = O (n^-'Y 



Decomposing in the eigenvector basis of m„ gives, with the notation of the proof of Lemma 

[HSl 



w 



M=^n. + E(^i"^^"'^'^'^")^'^i"Vr- 



fc=l 
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Averaging with respect to j G {0, . . . , n — 1} gives: 

n 
k=l 

where 

\\vn\\=0{ni-'y 

Hence, 

(n \ 1/2 

We then estimate the right-hand side of this equation, using the following lemma: 
Lemma 11.6. There exists a universal constant c > 0, such that for all X & V , n > 1, 

|M„(A)|2<l-c((n|A-l|)Al)2. 

Proof. If |A — 1| > 3/n, we have 

h _ \"|2 

\Mn{X)\'= ' ' <4/9, 

n^|l — A|"^ 

and the inequality is true for c = 5/9. If |A — 1| < 3/n, we can write A = e*^ where 

|A - 1| < 1^1 < 7r|A - l|/2 < 37i/2n. 
Then, 

^ n—l 

l-|M„(A)p = - 5^ (n-|j|)(l-A^-) 

j=-n+l 



= ^2 T. in-\Ml-cos{je)), 

j=-n+l 

where |j0| < 37r/2, which implies that 1 — cos{j6) > j^^^. Hence, 

n-l 



j=~n+l 

>n^e^>\X-l\^n^. 



n 



Let us now estimate (125|) . first by assuming that x is not one of the values of yk for 
k E 7j. Then, there exists k E Z, Xi,X2 ^^, such that 

Vk < Xi < X < X2 < Vk+i- 
Hence, for n large enough, 

—^ <Xi<X<X2< -1^, 
Zn Zn 

and there are no eigenangles of Un between 2TTXi/n and 211x2/ f^- One deduces that for 
all fc G {1, . . . ,n}, 

and by Lemma lll.6[ there exists d < 1 such that for all n large enough, and all k G 
{l,...,n}, 

|M„(Ai"^e-''"'^/")| <ci. 
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^M||<c/(X^l4"^n +0(n2-^) = d\\w[n]\\ + O (n^-^) 



Therefore, (1251) gives, for n large enough: 
|w|n]|| < d 

\k=l 

and, since 1 — (i > is independent of n, 

IkMII ^ (1 - ^)lkMII = o (n^-^) . 

By Lemma m. 5 [ w is identically zero, which implies that {V°')aeR has no eigenvectors for 
x^ {yk,keZ}. 

Now, let us assume that x = Z/fc for A; > 1 (the case A; < is similar) and let W7„ be the 
projection of w[n\ on the orthogonal of /^" . We have 

W[n\ =Wn + Vk fk ' 

and then 

where the two terms in the last sum are orthogonal vectors. Hence, 

sup I \u^{wn) - e2*">^^'/"w„| I < sup I \uiw[n] - e^^'^/^wM 1 1 < n^-^. 

0<j<n-l 0<j<n~l 

Now, decomposing Wn in the eigenvector basis of m„, and performing the same computa- 
tion as for w[n\, we obtain an estimate which is similar to ( 125|) : 

(26) \\w^\\=( J2 \MniX^^e-'-^/n\'\v^M +o(n^^' 

\l<fc'<n,A:Yfc / 

Notice that the term k' = k is not in the sum, since w„ is obtained from w[n] by removing 
the component proportional to /^ . Now, it is easy to check that for k' E {I, . . . n}\{k}, 
one has the estimate 

(note that this estimate is not true for k' = k, since nO]^ /2'k tends to x = Vk when n 
goes to infinity). Hence, for n large enough, |M„(A^? e~^*'^'^/")| is uniformly bounded by 
a quantity which strictly smaller than 1. As above for w[n], this implies the estimate: 

\\wn\\ = O (n-2~'y 
If we set Kn := rj^ /D^ , we deduce 

\\w[n] -Kngt'^W =0 (n^^'^J , 
and then 

(27) \\w[n] - Kngk[n]\\ = O (^{1 + \Kn\)n^~^^ , 
smcG 

Now, for any integer m such that n < m < 2n, we have 

\\w[m] - K.^gk[m]\\ < (1 + \Kra\)m^'^ = O (^(1 + \Km\)n^'^J , 
and taking the n first components, we obtain: 

(28) \\w[n] - K^gk[n]\\ = O f (1 + \Km\)n^-^ 



"^m I I "^n I _ I "^m "^n | _ r) V "i" I "■ I ' I "m \ ) j 
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Comparing ([27]) and (12S|) gives the following: 

|«m - K^lllsffcNII =0 (J,l + \Kn\ + \Km\)n^'^ 

and then 

(29) \Km - Kn\= O ((1 + |/t„| + |/€m|)^"^) , 

since ||5'a:[^]|| is equivalent to a strictly positive constant times ^/n. In particular, for n 
large enough and n < m < 2n, 

which implies 

|/tm| < 1 + 3|k„| = 0(1 + \Kn\), 

and (129|) can be replaced by 

(30) \K,m - K-nl = O {{1 + \Kn\)n'^) . 

Hence, 

\l^m\ = \fin\ + O ((1 + \Kn\)n'^) , 

sup {l + \^m\) = {l + K\){l + 0{n-')), 

n<m<2n 

and 

where Sg denotes the supremum of 1 + |fi;m| for 2'' < m < 2''"'"^. We deduce that the se- 
quence {Sq)q>o, and then the sequence {nn)n>i, are bounded. The estimate f l30|) becomes: 

for n < m < 2n. If for g > 1, 2'^n < m < 2'^^^n, we also get: 

g-l q 

\l^m - l^n\ < \l^m " l^2in\ + ^ |/«2'-+in " '«2'-n| < ^(2''n)"'' = 0{n'^). 

r=0 r=0 

Hence (k„)„>i is a Cauchy sequence, and if k denotes its limit, one has 

\k — Kn\ = 0{rr ). 
Using fl27|) . we deduce: 

\\{w - ngk)[n]\\ < \\w\n] -K„fi'yt[ra]|| + |fi: - /«„|||5'feNII 

< (1 + \Kn\)n-^-^ + n-^^ = O (n"^-^^ . 

Now, w — KQk is a sequence in J-": by Lemma [11 .S] w — Kg^ = 0, which implies that w is 
proportional to gk- □ 
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